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SYi-JOFSlS 

DESIGN OF A 120 MWt RE/aCTOR FOR A SURFACE WARSHIP ’ 

by 

LT. ARUN J. RAO, IN 

Never before has the strategy of maritime warfare 
diversifiea so much as with the ingress of the nuclear reactor 
onto warships. Subsurface or surface, the nuclear propelled 
warships are far superior both tactically and in performance over 
their conventional counterparts. With the fast dwindling number 
of warships due to their esclating costs, the possession of ! 

nuclear-powered warships for maintaining 'the same effectiveness 
and punch is becoming increasingly inescapable. 

Reactor design is a multi-dimensional complex problem- 
embracing many fields of engineering. The warship reactor design ! 
is more intricate than the land-based ones at it involves 
certain additional specialized problems like ship-motion effects, 
manoeuvring, vibrations, etc. The physical parameters of the i 

core is restricted mostly by thermal considerations. Hence j 

1 

the Thermal hydraulic design requires careful attention to 

1 

prevent core melt-down or coolant bulk-boiling. The Reactor 
physics design involves evolution of core enrichment and certain 
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other parameters to result in optimal design. Knowledge of 
the neutron flux distribution is also required. Here, two-group 
theory for heterogeneous core lattices has been used. The 
Shielding design involves meticulous calculations to produce 
the smallest possible reactor shield, both in size and weight 
which is of paramount concern in a warship design. 

Control of a marine reactor is different from 
the land-based ones. Due to stringent manoeuvring requirements. 
Power excxarsions from 20% to full-load require automatic 
control with a high power increment /decrement rate of 3 °^ . 

The system has to be prompt acting, stable at all times 
of operation with no unallowable overshoots. Finally, when 
the reactor goes on the ship, certain other problems like 
ship-motion effects, vibrations, etc., come in. 

Due to the immense calculations involved, any 
modern reactor design process is impossible without the aid 
of a digital computer. Here, 3 computer CODES have been 
developed, one each for thermal hydraulic. Reactor physics 
and Shielding designs. The CODES & their results have been 
presented separately. 



Tlie different design processes discusseu till 
now are inherently interactive among themselves and cannot 
be consiaered isolatedly. This has been highlighted in the early 
part of the thesis. It has, however, been the earnest t.nd.eu-Vv...ur 
of the autnor to ensure an overall non-conflicting design. 
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INTRODUCTION 


In the modern era of fast depleting oil^, coal 
and other forms of fossil fuels, nuclear energy takes a 
leading seat as the future global energy source. The nuclear 
reacfor is virtually an interminable source of energy. Tapping 
of this energy is an engineering feat by itself. Though 
marginally more expensive., it is only a matter of time in this 
energy-crisised. world. before the reactor, establishes itself 
as an inevitable source of energy. All this, and it is a mere 
42 years old 


Today the cost of a nuclear reactor is 
prohibitive. Why then nuclear propulsion? The answer is 
straightforward - No repeated- refuelling for the ship;, which 
opens an. entirely new era at sea,’ Where economic considar.ation 
is counterbalanced by military requirement., nuclear propulsion, 
for Naval ships become essential. For the merchant ships, 
however, this is not so. as money is the key factor for 

trade. Consequently, we have a, handful of Nuclear merchant 

whereas the number^ of navrrt c!v>-tr.,o .-v. , . 

propulsion) exceed three hundred. ^ world (with nucha? 

SCOPE OF THETHESia 


A compact low weight system is of ...paramount 
importance in plants designed for. ship propulsion. Marine 
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reactor design by Itself/ is a very vast topic merging land- 
based reactor design with certain specialized aspects like . 
marine environment/ light weight, limited space/ etc,^ Our 
goal is to design a reactor system for a frigate/destroyer 
class warship* A drawing board design involves a vasf. amount 
of dat49L,j tifne and generation of nximerous highly essential 
computer codes. This is definitely not being attempted. In 
the following pages, certain basic aspects of marine reactors 
have been discussed and a skeletal design attempted to 
highlight the varied & complex problem of marine reactor 
design. It is beyond the scope of this thesis to cover all 
topics envolving the design, be the skeletal stage itself.. 

All important aspects, have,, however, been covered. 

LITERATURE REVIEW 

Most of the information regarding Naval reactors 
io classified. Though information on reactoraoJ, merchant 
ships are easily available, they ..don't, satisfy the .stringent 
measures of manoeij^f ability, .high power-to-weight tatio. . 
and the limited space available in a wajrship where, economic 
considerations fall into the backgroun4»ik Yef/ a fair idea 
on problems enco,unter.ed in marine reactor design is avaiiat*le 
from them,. International symposims on Nuclear ships have , , 
been he.ld in (1960), n;an*>urg (1971) .and 

Madrid (1972). 
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A subcommittee on military applications of the 

U.s, has presented the tactical and effective use of Nuclear 

Propulsion for Naval Warships^^^ It argues that. though 

economically dear, when faced with the threat of submarines 

armed with either torpedoes or cruise missiles* or bombers 

or stirface ships armed with cruise missies / the all-nuclear 

task force (consisting of an aircraft carrier and a few 

destroyer s/frigates/ all nuclear— powered) is less vulnerable 

to attack than a task force of conventional ships. Ihe 

committee concludes# nuclear frigates for independent missions 

are far superior in performance than their conventional 

counterparts, John R Lamarsh^^^^ states that submarine reactors# 

^•35 ^ 

being very small and compact# use fuel enriched in u to .. 
over 90 w% A.N. Petrosy.'ant-s -e^^tj^uidating problems of 
nuclear science €£ technology# discusses Russian, ships with 
nuclear engines - 'Lenin' & 'Arktlka' (both ice-breakers) | 

j 

& their vastly superior performances over diesel icebreakers. : 

I 

I 

Extensive study & analysis of the world’s first i 

(5) . . ^ 

nuclear-propelled merchant ship 'Savannah # in. its sea [ 

I 

environment# has thrown light over the problems faced by. 

Marine reactors. Analysis of its core performance, has. .been, done 

by WJ. Gallagher and The Japanese Nuclear ship 'Mutsu* .. 

( 6 ) 

has also been discussed in detail by S, Sasaki Off-late# 

L. Chinaglia^^^ of Fiat# Italy has presented a new design of .. 

. 

ft 
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a PWR ,for ship propulsion. Ihe inherent design interactions 
among various design aspects like thermal hydraulic design/ 
nuclear design, shielding design, etc. has be.en well 
illiLStrsted:-. by ,RA Knief^^^^. 

El-Vakil deals with PWR thermal hydraulic 

design. Weisman ^ describes at length the various hot-^pot. 

factors concerning the thermal design. Regarding .reactor 

( 22 ) 

physics analysis, Meem , with his 2-gp theory, gives 

sample calculations for .reactor designers, that has been 

( 21 ) 

extensively referred. The Nuclear Engineering hancibook . . 

'^by Etherington has been used for related figures & 

references as and wl^iien mentic.»ned. IAEA reports Vienna, 

(27 ) 

1964 , highlight the reactor shielding design .and probl^. 

B, Chinaglia^^®^ of Fiat, Italy, has also worked on shielding 

studies Sc control of ship-based gaseous activity release, 

(29) 

J. Brown & A. Jackson , have carried out a study on 
shielding of Marine reactor installations. S.J, Ihunell^^*^^ 
has presented a comparative study on neutron fluxes in a • 1 

metal-watep shield in merchant ship reactors. 

i 

A M,arine Reactor control & in-core instrumentation j 

system design has been done by F„ Basile, R. Merino Sc Selviaoa. ^ 

! 

f 

wherein they have pointed out the diffetences between land-n 
based reactor, & marine reactor control. The control scheme i 

(2) i 

used in the thesis is one evolved by Chinaglia . of Italy 4 
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A detailed treatment of controllers has been illustrated 

by Schulte-: « this has been used, however, 

® study of non-linear control system. References ^for 
other aspects on reactor control have been Lewins^^*^^ and- 
Glasstone & Sesonpke^.^^^ and D.L, Hetrick^^^^. 

Specialized marine problems like ship-motion 
effects on steam-generator performance, collision barrier. 

£c grounding protection barrier for- ship s^re actors, and ^ 

safety of crew have been studied, respectively, by V, Jacob 
and G, Selvegg^^^^^^ W, Jager & H, lettnin and L. Chinaglia 

& P. Rucci^^^^ , 

1.3- THESIS ABSTRACT 

The thesis starts with a brief history of marine 
reactors. The tactical aspect of nuclear— powered vessels 
over their conventional counterparts is discussed- this may - 
m^e the difference between winning or losing a battle which ... 
is vital for any fighting force. The salient features of . 
ship-board reactors over the land-based ones.g^^ discussed. 

The modern trends o.f marine reactors ilso elucidated. . . 

Chapter .3 highli^ts the. design process, Ihe . 
different aspects of reactor design like thermal hydraulic 
design, nuclear ■ design, shielding design, etc., are inherently, 
interactive among themselves and hence cannot be independently 
considered; this is highlighted, A reactor .core .configuratioj:|i, 
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is then evolved and the systems used mentioned. 

Core thermal hydraulic design is developed in 
chapter 4, A linear power rating of 11.0 KW/m is used. 

Since the power generated by a reactor is essentially resti Ict^^ 
by the effective heat removal, this portion of the design 
is carefully considered as we desire the most compact and 
light-weight reactor for ship-board installation. A computer 
CODE, so vital- for any thermal analysis, has been developed 
and used to arrive at optimum reactor Parameters. With this 
knowledge, we next move over to the reactor physics design 
(chapter 5). Two-group theory is used and the core studied 
for criticality using 2 regions (water as reflector). A 
computer CODE has been evolved which computes the effective 
multiplication constant (K^) , the fast and thermal flust.es 
in both the regions, as well as the various 'poison ' (unwant-ed 
neutron absorbents) activities. in the core. The amount of 
control poison required for reactor control is also evolved. 

The control design hence would chiefly depend on the parameters 
obtained in this chapter. Care has been taken in the COpE 
that temperature and, pressure dependencies of all parameters 
ar e<'j:|;paj&idered . 

We are aware that the overall reactor weight 
and size has to be the minimum possible for a marine reactor. 
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BulSt of this weighty, however < is taken up by the biological 
shields’. Consequently^, reactor Yielding assumes considerable 
importance and is analyzed at length in chapter-6. Biological 
shielding is mainly required for the fast neutron and gamma 
fluxes emerging from the shield materials due, to thermal 
neutron capture, Ihe biological dose rate from the secondary 
shielding is to be restricted to 0.75 mrem/hr. This has been 
achieved after the formulation of a shielding CODE* In 
addition an approximate weight calculation of the reactor 
(with shield materials) is demonstrated. 

The, control of a marine reactor is greatly 
different from a land- based one. This is because of high 
manoeuvring demands and a requirement of controlling the 
reactor automatically .from full- load to part- load (20%).,, 
Chapter-7 elaborates these aspects and evolves a reactor control 
system. Careful stability analysis is also done from full-load 
to part-load using a two temperature feedback for the 
reactor. Certain other aspects of stability are also attempted 
after the calculation of the t«3^3er,ature coefficients of 
reactivity as applicable. A reactivity balance of the core 
using the CODE of chapter 5 is displayed. The control policy,.- ■ 
used ie elucidated subsequently. 
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The integrated aspect of the reactor with 
the ship increases the design complexity further. Its effects 
are of practical significance to the operator on-board, 

Chapter-8 embellishes this integrated view, Ihe ships three 
dimensional motion at sea has its effects on the steam 
generator performance. Collision barrier and grounding protection 
structxire is required for the reactor compartment. 

Chapter-9 concludes the present work 
illustrating the design that is ultimately arrived at. Other 
aspects of marine reactor design beyond the scope of this 
thesis are mentioned and suggest iv-ns for further study 
given. All references made in the current task are listed 
in the bibliography. 


* * * * * 
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In 1942^ Fermi & his co~workers made the world *s 
first reactor critical at the University of Chicago, However/ 
not until 1957 was a propulsion unit designed to operate on 
nuclear fuel. It was a small highly-enriched Pressurized 
Water Reactor (PWR), installed on board the Submarine USS 
NAUTllUS. As of today/ there are over 300 ships propelled 
by nuclear energy. Most of these, however, are warships 
where economy is not the deciding factor in the choice of 
propulsion system, 

2,1 BRIEF HISTORY 

2.1.1 NUCLEAR MERCHANT SHIPS 

The history of nuclear merchant vessels have never 
been bright. Nuclear power in merchant ships was used for 
the fir«t time on a 16,000 ton Russian ice-breaker 
in 1960, The decision was based not on economy but efficiency 
over the prevailing diesel electric system. In Arctic areas 
of heavy ice, nuclear ice-breakers with a 50% higher 

specific capacity*/ can operate most of the year; but the 
conventional ones can make way only for 3-summer months 
due to dangac of running out of fuel. The second ship was 
‘SAVANNAH', a U.S, Cargo-liner canmissioned in dec ‘62 
(no longer in active service). It's purpose was to test 
safety of nuclear plant aboardship & entrance of nuclear 

*SpecifiG capacity - power-to-displacement ratio. 
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vessels to civilian ports. 

The •OTTOHAHN* of ERG is a small scale experiro>sntal 
ship not planned for profitable operation (1968 vinf3‘ire). 

Her purpose is to furnish all types of information theori:. .eal, 
practical, economic & scientific with regards to nuclear 
propelled merchant ships at sea & refute arguments about their 
safety. The Japanese, on similar lines as JRG, launched a 
10,400 ton specialized nuclear cargo ship in 1970 (MUTSU). 
Russia, finding the nee! for a second nuclear ice-breaker 
commissioned ‘ARKTIKA* (displacement 23,460 tons) in 1975. 

The reason for this bleak growth of nuclear ships 

arc plain 

(i) Cost economy - Capital cost of nuclear ships are 60% 

higher than conventional. The fuel cost 
is also higher. Till such time fossil fuel 
becomes prohibitive enough to counter- 
balance the above two reasons, nuclear 
ships will not be' easily accepted. 

(ii) Uncertainties-of investment cost; the future development 

in this area is hardly predictable, 

(iii) Safety problems- Nuclear ships are not easily accepted 

by most of the civilian ports; psychological 
fear of public is mainly responsible for 


this 
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2.1,2 NUCLEAR WARSHIPS 

, A 

, , Nuclear Energy toox birth in warships in the 

form of a submarine reactor for 'NAUTILUS', a U.S. Submarine# on 
17 Jan '55, generated 15 #000 SHP# having 2 shafts and 

displacing 4040 tons. Initially- she steamed sutafnerged irom' Key 
West# Florida# to New London# a distance of 1397 sea miles?’, 

She was refuelled after 2 years and 62 #562 miles of steaming. 

Her success gave rise to a 'gold rush' and the US Navy ordered 
a pack of, submarine reactors from .Westinghouse Corporation. In 
1960# Russia started her sijbmarine fleet. 

'ENTERPRISE'# a 89#60C ) tons US Attack carrier,, 
was the first surface ship to have nuclear propulsion ,Cin 196i„),», . 
She has 8 PWRs on 4 shafts# giving her a power of 2^^ OOP , SHP. 
and a maximuni ..speed of 35 knots**. She has had her core char^geji, . . - ; 
thrice# with the present feed giving her a sustainence of 10 , 

to 13 years. In. 1962# ' BAINBRIDGE * # a U.S. Nuclear-Powered 
Guided Missile cruiser was the first cruiser to be nuc.lear- 


propelled (displacement 8580 tons). 

Changes in. =ship- reactor -technology have 
alsqi,l^n,,,inrmponated tp m^e ■ -ships - more versafiia- -andr:- ”-- — 
indapendsi^---!!*^ newer JUS., aircraft aarr, JJilMI^ (liajjgppjaeA.-. 


in ^,9,75) has pnly 2. reactors contrary to her' 



’ENTERPRISE* with^; this generates empower of .^0^#0pa ' Sf^ or 


4 st^fts.#., giving her ' 91.#40Q;^,tons a rnaximum speed-pf -.35- knot. 


She is the largest warship ©Y@r. , , 


* 1 ' se,e mile (or nautical mile) = 1.83 kms (2000 yariH^ ^ 

** 1 know =s i nautical mile/hour. 

\ - ' 
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built. She can steam an estimated 800,000 to 1 million miles 
before refuelling. Interestingly, the energy equivalent of 
her cores is 462 mega gallons ( niillion tons) of fuel 
oil- this is where nuclear energy scores over the conventional 
form. 

The largest submarine ever built, however, is 
a Russian one, a 'typhoon* class Ballastic Missile submarine 
(SSBN), launched in 1980. she displaces 30,000 tons dived 
with speed exceeding 24 knots and having 2 nuclear reactors, 
one per shaft, supplying a total of 120,000 SHP; she carries 
20 SLBMs ( Submarine Launched Ballistic Mxdsiles), Presently,, 
nuclear submarines fasE outnumber nuclear surface ships. 

However, many surface ships both cruisers and destroyers 
are increasingly becoming nuclear-propclled , mainly due to 
tactical reasons. 

Due to the prohibitive cost & superior technology 
for marine reactors, nuclear warships are expensive toys 
handled only by a few econcmically advanced countries. However 
with the Indian Ocean swarming with subsurface surface, 
nuclear & conventional warships, possession of nuclear warships 
to effectively safeguard our exclusive economic cone . & 
sea- frontiers, is becoming increasingly inescapable. 
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2.1.3 FUTURE TRENDS 

For merchant ships greater than 70,000 tons 

, (23) 

displacement, preliminary studies indicate nuclear Power 

to sc ,re over conventional power when operated for long 

durations* This might be used for constructing huge nuclear 

container ships in future when nuclear energy gets accepted 

in the merchant marine field. The smallest nuclear ship ever 

built is a US nuclear -powered ocean engineering & research 

vehicle 'NR-l* (1969), a sulanersible, built by USAEC, 

Knolls' Atomic Power laboratory. The vessel displaces only 

t I • 

400 tons, having dimensions of 136,4 x 12.4 x 14.6 with 
2 propellors & carrying 5 crew members and 2 scientists. 

Such vessels are very useful for conducting detailed 
undisturbed experimentation on the continental shelf over 
long durations. 

Now we hear the concept of NUFIS - Nuclear 

Floating Island, brought forward by Japanese scientists. 

Because of small area & dense population, the siting of 

Nuclear Power Plants in Japan is limited. Hence, the feasibility 
(9 ) 

of NUFIS has been studied and the project is about to 
commenu®. it is to be located 20 kms offshore at a sea-depth 
of 150 m. It is to have a 1200 MW^ PW Plant on a 140 m x 140 m 
floating platform with a displacement of 298,000 tons. The 
power plant &. the associated facilities are to be situated 


4 
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on it. Hence, the marine reactor technology will have to be 
applied to it. The technological feasibility of a floating 
centre for fuel fabrication, reprocessing, etc., is also 
being thougtit • 

2.2 TACTICAL ASPECTS OF NUCLEAR WARSHIPS 

It was stated that the elimination of repeated 
refuelling opens an entirely new era at sea. Let us explore 
the tactical boost this gives. 

Nuclear propulsion enables warships to traverse 
across seas at high speeds without peiising to refuel and 
without repleinishment. With the escalating cost of warships, 
many many times that of equivalent merchant ships, the 
number of warships in the fleet is bound to dwindle. Hence, 
'nuclear propulsion becomes prominent to maintain the same 
effectiveness with lesser warships. With the fast depleting 
oil sources^ in future, nuclear propulsion is going to be(pome 
a necessity. It gives an additional tactical flexibility ' and 
opens new courses of acti-,n to the operators of the Naval-., 
Forces which would otherwise be unavailable, A nuclear .> 
powered naval force will be able to achieve many more i 
productive results in any operation; such a force would) be, 
in addition, less vulnerable to the action of the enemy, W 

*' V 

/ 
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2.2.1 DRAWBACKS OF CONVENriONftL SHIPS 

Here, mobility is limited by the need for refuelling 
at sea. Enemy action, bad weather & other uncontrollable 
' forces directly impact upon refuelling ability which is not 
only time consuming but also requires the services of a 
considerable © 43 . tanker fleet/train. The Fuel oil carried by 
ships being inflammable, in addition, poses a fire .tte.o 2 ;ard 
onboard. A nuclear fleet, contrarily, is freed from the 
constraints of tankers & base support; hence it has increased 
speed, reliability, range, less fuel requirements & more 
payload . 

2.2.2 NUCLE/vR SUBMARINES 

In the modern warfield, no other weapon has proved 

to be such an imminent threat as the nuclear submarine::. As 

we see from Table 2.1, the super powers have maintained a 

peirity on the number of nuclear submarines. Admiral Moorer 

( 4 ) 

of US Navy has stated that he would rather have one nuclear 
powered - tactical cruise missile submarine (S^BN) than 5 of 
its conventional counterparts. Submarines are a particular 
threat because of their potential for approaching to within 
weapon ranges without being detected. 

When their batteries run down, diesel submarines 
have to periodically ‘snort?'. The quantity of fuel they carry 

*'Snort*“ ccxne up to near sea—surface level to to run diesel 
engines for charging batteries that are used for 
propulsion underwater. 



limit their edurance. The nuclear submarines^ ostensibly, ' 
never require to surface and can continuously move at high 
speeds cf 20 — 30 knots ( almost double of diesel submarines! )• 

The modern nuclear submarine being more quieter*, have a 
high probability of going undetected by enemy sonars. Hence,'' ^ 
they are much more effective in penetrating anti-submarine 
defences than diesel submarines. 

2.2.3 aURFJ^.CE SHIPS - THE ALL NUCLEAR TASK FORCE 

A naval task force consists of an aircraft 

carrier with her modern air-wing, and her escorts of 3 to 4 

destroyers/ frigates to provide additional protection against 

surface, sub-surface and air threats. The capability and 

punch of the force is restricted by the endurance of each 

of its ships. Hence, fcxr maximum effectiveness, all ships 

in it have to be nuclear-propelled . When faced with the 

threat of submarines either armed with torpedoes or cruise 
missiles, or bombers, or surface ships armed with cruise 
inifl^il^a, the all-nuclear task force is less vulnerable to 

attack than a task force of conventional ships which must 

show down to refuel in an extended encounter. It is therefore 

much more difficult for the submarine to get in a favourable 

position to attack the all-nuclear force. Further, the 

logistic support forces upon which the conventional striking 

forces are more dependent, are even more vulnerable to attack tha 

* ^quieter* - Vibrations of hull due to inner machinery is 

of lesser magnitude, making passive underwater 
detection more difficult. 



warship^. 


The tactical advantages of nuclear propulsion for 
ships of a nuclear task force can hence be summed up as 
follows 

They can move where required at top speeds without 

establishing logistic support before departure. 

Their virtually unlimited range enables them to use 

their superior speed to arrive at crisised situation 
much earlier than their conventional counterparts, 

When they do get there in a less vulnerable condition/ 

they have far greater staying power, not only in terms 
of their ability to stick on without fuel, but also, 
significantly, in terms of aviation fuel & aviation 
ordanance for their airwing. For example, a. 

‘NIMITZ* has 50% more aviation ordinance & twice as 
much jet fuel than equivalent conventional carriers, 

- It can continue combat operation without support 

until aircraft fuel, ammunition, food or other 
supplies are exhausted. 

It can disengage if required to retire at top speed 

to less vulnerable sources of aircraft fuel & ammunition 
replenishment. 



2.2.4 NUCLEAR PROPULSION FOR FRIGATES 


By virtue of possesing their capability for high 
sustained speeds and uninhibited endurance, nuclear frigates 
are uniquely suited for assignment to independent tasks with 
an effectiveness no non-nuclear ship can match. This can mean 
difference between winning or losing engagements. In independent 
operations, it matches the endurance of any nuclear sutmarine. 
Since the aim of the thesis is design of a reactor for frigate / 
destroyer, let us explore its tactical aspects a bit deeper. 

Nuclear frigates combine the essentially unlimited 
endurance of nuclear propulsion with the capabilities of the 
most modern anti-air and anti-submarine weapon systems. Their 
nuclear propulsion results in major improvements in offensive 
& defensive capabilities 

greater attack effectiveness due to increased mobility 

& freedom of independent action; ability to transit at 
high-speeds and thus arrive in attack position earlier 
& ready to conduct operations immediately without pausing 
to refuel; ability to be on attack station a higher 
percentage of time; the freedom to extend attack to 
larger areas, 

reduce vulnerability due to freedom from dependence upon 

replenishment in contoat areas; and ability to avoid 
submarine attack by transiting at high speed using 
evasive transit routes. 

significantly reduced nee-d for sea-based & land-based 

logistic support forces. 
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;- 4 ncareased propulsion plant reliability duer-to the 
higi^er standards required in the desigii <>5 con struct ion.., of 
nuclear propulsion plants. 

increased overall effectiveness of the entire fleet 

by reducing the total ntimber-.of coirbatant. ships„re.qvtJxe< 3 , 
to protect , oilers , 

2,2., 5 WORIO'S MUChEAR FliEET STRECJGTH 

Table 2.1 lists the nuclear strength of the world 
navies as of today, Ihough nuclear submar^es for outnvutiber 
nuclear surface shlps^ we can see that, the nunhej^of the 
latter has more than tripled in the last, 10 years and many 
are being bnilt» 
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,Mj^lNE R5ACTQRS 


L/'J^D BiwSED REACTORS 


The crew (& passengers) 
constitute a ‘captive* augiienco; 
they have no place of shelter 

in case of accidents. 

In such an event/ • 

people could be evacua- 
ted to safer areas. 

Most likely place of accident ! 

is in the least favourable 

location - the port areas. ; 

This problem is 

non-existont , , 

, (ii) DESIGN 


There are many more start-up 
and shutdowns as well as widely 
varying needs for power in man- 

csauvering. This results on a 

greater stress & strain on the 
nuclear machinery (to be taken 
care of in .design. ) 

The design need not 

be that stringent. 

1 

Due to space limitation the 
machinery of working plah-fes -^'9 
much more crowded than in 

stationary plants. 

No such limitations,^ 

. Due to weight & space limitat- 
ions'/ the shops / spare-parts 

and maintenance facilities are 
• 

kept to a minimum. 

- do - 

Design has to cater for specialist 
aspects on ships like vibrations 
(of propeller & machinery) / sea 
conditions (roll/ pitch/ yaw) and 
limited space and weight. 

Earthquake protection 
is required for reactors. 






25 




REACTORS 

11 . 

Irrespective of economic I 

Materials chosen should 


considerations, only the best . : 

be such as to minimize 


materials are to be chosen such j 

capital & maintenance 


as to achieve highest possible 

] 

cost while still being 


specific capacity. ! 

effective. 

12. 

Reactor control should be such 

aa.to make it inherently safe, i 

in the event of a tsvaiamd or 

sinking of the ship, 

i 

, No such restriction 

13. 

i 

Propulsion system design should 
cater for an emergency conventi- 
onally powered propulsion system! 

No such condition. 


compatable with the main machinery. 


in the event of a complete break- 
down. j 

1 


14. . 

j 

In warships, adequate shock j 

No such condition 


protection for reactors is required required to be 


against firing of onboard gund 
and missiles? also required .is 

(i) protection against shelling 

(ii) grounding & collision 
protective!! against external 

accidents. 

catered for. 

15 *, 

Only PWR * s are used as marine 

All types of reactors 


reactors. 

are possible on land, . 
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2.,3».2 CRITERIA FOR -WClgLj!^ M^SHIRS 

Thre.ats to crew health & safety do not arise 
out of normal operation of nuclear ships hut rather du,e 
to the possibility of internally and externally paused 
accidents to them. To minimize these# a variety of special 
safeguards in design# construction & operation are required. , 
These will be briefly elucidated. There are also# a variety 
of recommendations applicable to nuclear ships; these# 
however > won 't be discussed here but can be found in 
references (7) & (10), 

Considerations such as the following lead to, 
certain basic objectives that must be achieve^- in the design 
of a nuclear warship; 

- There mxist be contij^ual knowledge of the level of radio- 
activity in, all normal effluent streams.#. & of the Iqyels . . 
at various -locations around the ship, 

- Ihere mu.st be complete control, of .all normal effluent 
streams,. ,w,ith storage capacity sufficient to permit detailed 
release of effluents under unfavourable circiimstances,^ 

- Ihere mtust be- an extraordinary degree of conservatism/ 
redundancy & dependability in the design of the reactor, 

so that the likelihood of .accidents will be. redpoed, -tp t^ 
Ipweet practicable level, following -aspects ere importaiits- 
^.l^arine environment like roll# pitch#. yaWf 
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^ PrQpeller & machinery vibrations of the reactor 
donq^sartinsnt . 

Proper reactor containment design for safety. 

Collision & grounding barrier protection for 
reactor compartment. 

/ 

Protection of reactor against shocks impressed 
by weapons firing onboard. 

The ship should have highest possible specific 
capacity. 

- There must be devices & safeguards to minimize the hazardous 
consequences in case serious accidents should occur. 

- There must be an adequate conventional emergency propulsion 
system in case of complete main-engine breakdown. 

Thus the importance of nuclear propulsion for 
warships as a means towards strategic mobility is obvious 
to all who pause to consider the matter. It will remain as . , 
until a superior form of ocean propulsion for the weapon 
sys1::ems is produced. 


Vf ★ ilr ★ ★ ★ ★ yr 
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It is well known that only PWRs are used as mobile 
reactors. In this chapter, this has been justified. In 
addition, the design methodology, the fuel assembly design 
and the . steam-cycle used has also been elaborated. 

3.1 CH OICfi OF PEACTOR 

■“Si-iyiCOLF" , ^ second nuclear -powered U.S. Submarine 

was fitted w.ith a liquid metal (JSIak) reactor called the 
Submarine Intermediate Reactor (SIR) in 1956. Dxxring dockside 
trials, steam-leaks deveoped due to Nak alloy coolant entering 
the superheater steam-piping. After repairs' she was put to sea. 
However, due to reciirrence of the ajDove leaks in Dec '58 , 
the reactor had to be replaced with a PWR. 

Ihe above experience was an eye-opener for the choice 
of mobile reactors. Alternately, if BWEJ were. used, they would 
give inherent system instability due to the reactor water 
level constantly varying with the ship^motion effects; in 
addition, for direct cycle the turbines require shielding, 
hiking up the gr^ss.. machinery weight which is a very undesirable 
phenomenon. 

Here, it is not possible to compare the various 
types of reactors due to space limitations. However, the 
salient features of a PWR over other reactors have been listed 
below:- 
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Advatotagaa 


, - Moderator and coolant same. 

“ Water easily available ana cheap* 

“ Weter poses no serious long-term radioactive hazard. 

- A primary h^at exchanger leak is not hazardous to 
personnel as. primary loop is shielded. „ 

- Reactor has a high negative temper atvire coefficient 
making it safe & stable to operate. 

_ For moderate enrichment and closely— packed lattices^ 
burnup is the highest. 

- Being pressurised, the reactor is compact and 
size a minimum due to high power densities. 

- technology is very well known as it is. the 
rat.. St common tip's of naval react<jr tcday. 

Pis adv.ant aues 

- Hi^ primary loop pres stir e required due to the low 
l^aJ-ing pQ,int of water; so pressure vessel should 
be very strong. 

-Poor thermodynamic ef f iciency-20-25% for mobile 
reactors. 

- As, ordinary (light) water has a large neutron 
absorption cross-section neutron economy is poor. 

- Reactor should be shut-down for reloading.; 

( i.e, reloading not possible at sea). 
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Despite its npderate disadvantages, the power-to- 
weight ratio, of a PWR is the highest; this is a highly 
desirable, , factor where compactness and light-weisbt is .a chief 
requirement* All marine reactors operating today are hence 
only PWFs . . 

3.2 THE DESIGN SCHME 

3,2.1 THE PROCESS 

This may be divided into the following categories:- 

(i) Core thermal hydraulics design 

(ii) Core Nuclear design 

(iii) Shielding Design 

(iv) Reactor control and safety 

(v) Materials Design 

Since each of the categories is complex by itself, 
it 4fS - being treated separately in chapters that follow. 

The last topic (v) , however, has not been included aepfiarateJ-y^ 
but sketchily covered throughout. Being a reactor for defence 
purpose, no study on cost -economics has been done. . 

3.2,2. INTER-^TIONS 

The integrated reactor design requires substantial 
interaction among often conflicting design principles of 
(i) to (iv) abovej safety considerations like biological dose , 



DESIGN SPECIFICATIONS 
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lin^tatlori/ prevention of core melt— down/ etc, provide very 
stringent constraint&4 

Fig 3,1 gives a pictorial view of these interactions 
in the present design. Nuclear/ thermal-hydraulic and 
material design components result in composition/ geometr;^' 
and operational specifications that characterize the marine 
.reactor system. Nuclear design begins with preliminary 
cross-sections data. For the calculational modules, flux- 
averaged cross-sections with appropriate modifications 
to the working conditions have been used. The fliox, power 
distribution and biornup contribute to a final specification 
of enrichment. The thermal hydraulic design interacts with 
power distribution to develop a geometry and its related 
temperature and flow distributions, limited by material 
specifications. The nuclear design also yields the reactivity 
requirements required for the control design, 

Ihe control system design for marine reactors is 

intricate. Whereas a power variation of only 10% A@»ijal<5'2>wed, 
in land-based ]?WR ' syar'iatiun from 100% down to ' 

20% is required to be automatically controlled that too in 
much shorter times for marine reactors. This calls for a 
non-convent ion al control system design with careful stability 
analysis. This is dealt with in chapter-7, 

3.3 FUEL ASSEMBLY DESIGN j The first step in core design 
is to evolve the type of fuel assembly required 




Fig. 3. 2 

The fuel assembly 
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The basic features of the assembly selected is shown in 
Fig 3,2 (Ref. (7)), A 17 x. 17 square pitch layout has been 
chosen with 282 cylinderical fuel rods per assembly. This 
conf igiiration represents recent designs that enhance linear 
heat-rate limits ( Ref, (13)), The taxonomy of the fuel 
assemblies is as shown in Fig, 3.3 . This symmetrical 
arrangement employs an open lattice that permits some flow 
mixing between adjucent units which is a desirable phenomenon,. 
The diameter of the core baffle is 7.071 times the width 
of a fuel assembly, 

A constant fuel enrichment has been chosen 
throughout the core. For a linear power rating chosen as 
11.00 KW/m (see chapter 4) and 120 MWt power output/ the 
number of fuel assemblies required are 29 (8178 rods). This 
gives a core height of 133,39 cms. The core-radius is 
dependent on the fuel rod-pitch (which will be determined 
later) , 

6 ' 

There are/control-roci poison fingers per assembly- 

each finger being of the same length the aose, .It cpntains 
B^C ( Boron carbide) enclosed in stainless steel tubes. Their . 
concentration and size are to be determined. The reactivity 
control is provided by 29 Rod Cluster jCcoatroi (RCC) assemblies 
(1 rod-cluster per assembly). They have been constituted in 
3 banks A, B and C as shown in Fig, 3,3. No btornable poison 
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or soluble poison scheme has been included. Each assembly/ 
in addition, has been provided with a central instrument thiml;!^ 
space meant for a movable in core flux meastirement detector 

3.4 STEAM CYCLE 

Ihe steam cycle used is illustrated in Fig, 3. 4, 

The 2- loop system is in vogue in .almost all working PWR's. 

The; primary system contains water ^t 2000 psi (14 MPa ). 

Coolant pump regulates flow in the loop and pressure vessel. 

The presurizer maintains the system pressrore within a specified 
range. ^,^2 primary coolant loops supplying to 2 u -tube ste.am 
generators of 60 MWt capacity each are chosen, 

■ i The secondary loop generates steam at 500 psi 

(3.5 MPa) at 100% load at a qpaality of 0.9975, Regenerative. 

been 

cycles for the turbines have^useci -to improve efficiency, 

Fxnrther details about the steam cycle have not. been attempted. 

3.5- INFERENCE 

We have now determined few of the reactor core 
parameters and the system we have to work on in the. succeeding 
chapters. We also know about system interactions. V^Jhile , 
considering the various design processes in the subsequ.en|; 
chapters, these interactions have to be borne .in mind for 
an overall non-conflicting design. This has been the earnest 
endeavour of the author. 
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4,1 INTRODUCTION 

A nuclear reactor is capable of operating at very 
high power densities. However, the core must be operated at such 
a power level that with the best available heat removal systems, 
the cladding and fuel t emper at ur e s anywhere in the core do not 
exceed safe limits. A fuel element damage, otherwise, will result 
either in the transfusion of the highly radio-active fission- 
products into the coolant, or in core- fuel meltdown. Both these 
effects are highly ■uiri'v^esirable Hence, the core design is 
dictated more by the thermal limitations rather than any nuclear 
cons ider at ion . 

GOVERNING EQUATION 

The rate of heat removal, (Q) is given by the simple 
relation (p:369, (17)) 

Q = h * A ^ 

ov^all''''''''''^Area of Tempr.diff. between 

ht. transfer heat transfer fuel & coolant, 
coefficient* 

'h ' is governed by conductive heat transfer through fuel elemerits 
and cladding, and convict iv© heat transfer through coolant. As 
the fuel is UO^ , its thermal conductivity is fixed. So fuel-; 
cladding should have least resistance to heat conduction, in 
addition to low neutron absorption cross-section, workability). 
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ability to withstand high temperatures and withstand the ciorr<o^sive 
action of the coolant fluid; some of these considerations are 
conflicting; hence a compromise between high thermal conductivity 
and high temperature operation necessitated the use of stainless 
steel (^^-347) as the cladding material. Resistance to heat 
conduction is reduced by keeping the cladding thickness to the 
minimum. Heat transfer by convection to water depends on the 
properties of water at the operating conditions# and coolant 
velocity. 

"A" is max-irnized by having as many fuel rods 
as possible. Since a large quantity of heat is generated in the 
compact core# the maximum possible heat transfer area should 
be squeezed into it. 

"/Nt"# finally# should be high; high rates of heat 
release requires low coolant temperatures; this conflicts with 
basic thermodynamic coJitoiderations as a good thermal efficiency 
requires high coolant temperatures. Hence# to, kepp coolant 
temperatures as high as possible and yet retain a high At in the 
core, fuel elements have to operate at the highest possible 
temperatures which are limited only by metallurgical & burn-up 
considerations. 

Hence# it is seen that reactor core design necessital 


a compromise between often conflicting nuclear# metallxirgical. 
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structural^ thermodynamic and heat -transfer requirements. To 
ensure that the core is safe under all operating conditions 
and transients in our main motto. 

4.2 DESIGN CHARACTERISTICS 
4.2.1 SALIENT FEATURES 

The main problems encountered in thermal hydraulic 

design are 

(i) assure inlet & outlet core flow distribution compatable with 
power distribution foreseen during the core lifetime. 

(ii) assure cooling of the control rod pins and hydraulic dashpot 
action (for slowing down control rod velocity at the end of 
its insertion) during all operating and accident conditions. 

(iii) DNB (Departiire from Nucleate Boiling) should not occxir 
under all operating conditions & tramsients ( i.e. Heat 
flux anywhere, to be kept below CHF (‘^ritifal Heat Flux)^ 

(iv) Bulk boiling of the coolant should hot be allowed to occiur 
in any channel. 

(v) The fuel centre-line temperature should not exceed 2500°C 
under any circumstancs (fuel m.p. 2850 °C ) 

(vi) Hot-channel factors should be judiciously chosen; with 

the heat-flux increased by these factors, points (iii),(iv) 
and (v) above should be satisfied. 
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4.2.2 CORE FLUX DISTRIBUTION 

The basic core flux distribution has been taken 
from the "Reactor 'Physics Design" section (Ch,5). Radial flux 
distribution has been assumed as a truncated Bessel's function; 
the axial distributioUj truncated cosine. Due to the high inbuiJ't^ 
reactivity of the core,, control rods remain partially inserted in 
the core. At the BOC ( Beginning of Life/ of Core) tbe flux peaks 
as shown in Fig, (4.1)/ (A). As the rods are gradually retracted 
due to burnup/ the flux assumes a more sinussidal distribution. 
Finally at the EOL (End of Life)/ the distribution assumes profile 
(C). For calculational purposes/ (B) is assumed as Sin(Z) 
and (A) and (C) as Z-Sin(Z) and ( ) . sin ( •tt-Z ) respectively. 

For comparitive purposes/ the normalized flux distributions have 
been shown as a computer printout. ( See p.27 of computer printouts) 

4.2.3 HOT CH/vNNEL F/xCTORS 

In the thermal design of a marine reactor core one 
is concerned with the most severe conditions in the core. If 
the design is adequate for these conditions/ it is certain to 
be satisfactory elsewhere. Hence the 'hot-channel ' is that coolant 
channel where heat-flxox and enthalpy rise are at a maximum. 

(a) Nuclear hot channel factor (Fq ) 

or Nuclear heat flux factor is defined as: 



44 


Estimated maximum heat flux in core (based on nominal 
“ fuel enrichment & dimensions), 

Mean Heat flux in core 

* F^ where, subscripts 'Z' & 'R' stand for axial 
& radial flux factors. This is the peaking of the fltix in the 
axial Si radial directions given as "leaking Factcrs'Mn the CODE. 

F^ < 2.32 ; F^ cC 1.57 for Sin(z) fl^ax to Fg^' 1.82 

for Z-Sin( z ) and ( 7t-z) gin.. jp^-Z) fluxes. Tnis causes 

F^ .^4. 222.4 ( We have obtained F^ = 3.4411 for Z‘Sin(Z) flux,) 

This factor has already been included in the CODE. 

E 

(b) Engineering hot-channel factor (F^ ) . 

The maximum hdat flux may increase under the ^ 
following engineering uncer taint ies;- 

(i) local fuel density being greater than nominal (factor 

E 

concerned - F,^ ) . 

E 

(ii) Fuel dimensions being in excess of nominal (F^ ) 

— 

(iii) Local fuel enrichments being greater than nominal (F )• 

e 

F^ — 

For a PWR this gives a value of 1.08 to 1,10 
(Ref (26), pp. 170-177). But this approach is unnecessarily 
conservative and the effect of fabrication tolerances is evaluated 
statistically giving 

Fq = 1.04 . (^Ref.(20)) 



45 


E 

(c) Engineering hot-channel factor for enthalpy rise )• 

This "accounts for engineering uncertainties in 

the maximum enthalpy rise in the concerned hot-channel. 

E 

depends on the follovjing f actors ;- 

(i) Fabrication tolerance subfactor, accounts for fuel j 

enrichment variation, fuel weight variation & coolant flow , 
variation; it is evaluated statistically. 

(ii) Core Flow Sub-factor (E- ). accounts for the hot-channet 

Ah , Flow 

receiving less than its expected share of coolant flow, 

E 

(iii) Core mixing sub- factor (F, , . ) , accounts for flow from 

uwi,.niix/ 

surrounding cooler channels mixing with flow in the hottest 
channel,' This is favourable countering the deleterious 
effect of the other subfactors. 

Typical values of above for a PWR are:- 

,E ^ 


* p'^ * p* 

■^yO-h Ah, Fab Ah, Flow Ah, mix 


1.08 


(d) Overall hot— channel Factor (F) 

E 

The engineering (F ) uncertainty factor would be 

F® = Fq * ( 1.167 for PWR ■) 

However, in the Italian ship Design (Eef(2) ,p. 598) 1.22 has 

E 

been used instead of the above figure, giving F =l-,'269; this 
has been used in this analysis. This means if the heat flux 
is increased by this amount in the CODE, the core design siaou-id 



should sustain itself. (F^ has, as seated earlier, already bee 
included in the C' i: ). 

The overall factor F, would hence ba; 

N E 

F = Fq * F = 4.357 

This implies the core should sustain a heat fl\ix over its a'^-'.r 
value, of 4.367 times everyv/here in the core. 

4.2.4 THE THERMAL HYDRAULIC C ODE 

The thermal design of a core is virtually 

impossible without the aid of a computer. In the computer 

care has been exercised to entertair all thermal effects end 

variation of material properties with temperature. It lists 

output in a 3-dimensional manner giving axial variations of 

temperatures & pressure for the rods at the core centre, an'": 

then at discrete steps upto the core periphere. Ihe analysis 

of the rod is for;- 

(i) Coolant integrity 
(it) Cladding integrity. 

(ii'O Coolant exit conditions and pressure drop. 

In addition for:- 

(iv) Radial fhrx variation in the core 

(v) P’omping requirements for primary system. 

(vi) Determination of core diameter (for 120 M4th power), 
and total number of rods. 

(vii) Temperature feedback for reactor neutronics calculation s. 
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The importance of the CODE is felt when one has 
to do repeated analysis of a core with small variations in 
input/ to arrive at an optimum solution. In addition/ cores 
of different reactors can easily be analysed since material 
properties are fixed. The CODE developed is versatile enough 
to be used either for stainless steel or Zirconium cladding, 
and either a square- arrayed ■ core or a hexagonal arrayed- ’ core. 

4.3 THERMAL AMD HYDRAULIC DESIGN 

The fuel rod has been axially divided into 50 equal 
segments. The properties of materials in one segment are considered 
constant ana so also temperatures^, axially. The basic theory , 
governing the temperature distribution/ heat transfer rates, ;j 
pressure drops & energy balance is presented below, and done 
for each of the 50 segments, successively. Temperature distribution 
in fuel Sc cladding are calculated using one-dimenS ional radial ^ 
heat conduction. -The pressxire drops & heat transfer coefficients . 
are calculated using empirical relations for flow channels. 

In this core design, the square 
arrary has been used. For the 
fuel rods of diameter DCL an<^ 
pitch s / the flow ape a (A) p^r 
rod is given by | 

A =5 s *DCL^ (shadea portion^ 

. . . .Eqn. .(4.1 ) \ ' 


4.3.1 FLOW CHANIJEL 



4 
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4.3.2 HEi*T GENERATION 

Fission heat generated is directly proportional to 
thermal flux (taking cross-sections as constant). The volumetric 
source strength 

* * * O 

^ “4 = 3.204 * 10”^^ *• “'g ' J2f (r,z) . J/cm -s 

* >i ^ 

Thus, q varies axially & radially as: 

q''(z,r) =q_ * Cos(jc z) * J ( 2.405 .r) Eqn (4.2) 

° R+d 

«. 

for 2=*0 to' H 

I I t 

r=0 to A-lt^vhere is the peak flux at centre of core, 

R= Core Radius 

d= distance from core face where 
flux - 0. 

Jin important quantity required to be known beforehand is heat-flu:, 
generated per unit length of rod (kw/m ) q^^ . This is an average 

quantity for all rods. After a review of the reactors of the 
world (Ref. (14)), a figure of 11.0 KW/m was arrived at. 

This factor is related to the average volumetric heat flxox thus- 

I 

/ti where, r£^ = fuel pellet radius 
The volumetric power density at the core centre is •- ■. l 


1 1 I 



■R 


★ ★ 

^2 ‘ 

,N 


I I I 

q 


whert; F^ & F 
before. 


av 

,N 
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eire hot-channel factors as explained 



N :s po z 
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since by virtue of definitio^v 






4 - fl ■ - ■ 

« ■! J J' 

V all 

o 

core 


Cos (JLz ). 
He 


( 2.40 5 

R+vi 


2 • rdr^lqe 


where V= core volume 
we arrive at the following relations 


2 

mR 


for a truncated Bessel ' s 
function flux-distribution^ 


“TT *— rrr / Sin ( "t „) , for truncated cosine 

*yrr 

”e 

1.8197 * ^ * H 


He i-hCosiA) + B Cos(B^+ Sin(A)-Sin(B) ) 
^ truncated ( t-z).sin (t-.z) function 


, for 


where y A = H Tt 
He 


and B= (' 


He-H Tr.t 
2He ^ 


Refer Appendix A3- 1 for the above inferences. 

4.3.3 temperature DISTRIBUTIONS 
(l) In fuel rod - 

The basic heat conduction equation with heat generation 
is (in cylinder ical coordinates ) 

, 1 dT 


2 

d T 


dr 


dr 


I I I 

Fu 


= 0 
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with boundary conditions as 

T=To at r=o (fuel centre-line temperature) 


dT 

dr 


= 0 at r=<o 


the solution of the above equation is 

I t I 


= To 


4K, 


Fu 


showing a parabolic temperature distribution with the maximum 
at the fuel rod centre. Hence, at the fuel outer radius (r^,^) , 

I t t 


Tp = To 


q 


4 K. 


'Fu 


.... ( Eqn 4.3 ) 


Fu 


During the steady state condition, all of the energy generated 
must be transferred out to the cladding - since changes with 

temperature, an iterative technique has been used to obtain the 
cladding thermal conductivity to within 1%. 


If K^(T) = a+bT where a and b are constants, 

then we find that the value of to be used is precisel 

To + T 

equal to K evaluated at the mean temperatvire F 


Fu 

i.e. in Eqn. 4.3, Kj^(T) = 


(Refer Appendix A4-2) 


(2) In gas gap 


The temperature drop across the gas gap is taken as 

( Eqn 4.4) 


T^ -T = 
f g 


r-, 

J 


2 h. 


gap 


ICl 2 

where gap = gap conductance = 0.6w'^^m -K 

(p,356, ref (12)) 
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h is taken as constant throughout the core, 
gap 


( 3 ) In the claddinc 


Technique used here is similar to that in the fuel rod. 


The governing equation is :■ 


1_ .±_ (r ^ = 0 ... Eqn.(4.4b) 

r dr dr 


with initial conditions as 


T = T at r=sf^ (outer radius Of gap ) 


^ 'I 


r = r. 


where, thermal conductivity of cladding. 


Eqn. 4. 4b gives a simple logarithmic solution 

Ml 2 ^ 2 

T = T.. - " fu ^ / r_, 

g — r— p + ln( Cl ) . 


(Eqn4.5) 


where r t outer cladding radius, 
ci 

Here again, is a linear function of temperature and hence 

an iterative technique has been used to obtain to within 1% 


occur acy, 


K ,(T) = K ( '^F'^ '^S ) 

cl d 2 


(4) In the coolant 


All heat energy generated by the fuel must be ultimately 
transferred to the coolant through the fluid layers near the 
cladding outer surface. Newton's law of cooling has been used 
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for this purpose 


/ 

heat transfer heat transfer '^heat transfer 
rate. coefficient area 



Tg is a mass-weighted average tempc-rature of the water in the 
flow channel. 

As, 


A 


s 


2 r^^ 


thuS/ 



» i 1 


Fu' 


2- h r 


ol 


(Eqn 4.6 ) 


The method of determining (h •) 'h: elucidated later in this 
section (4.3.5). 

Tg is to be known 


before we start temperature calculations. After that. 

To are determined using equations 4, 6, 4. 5, 4. 4 and 4.3 

successively. Method of calculating Tg has been stated later 
(4.3.6) . 

4.3.4 COOL;iI\lT PRESSURE DROP 

We shall determine only the pressure drop in the core 
here. There are 2 types of losses occuring. 

(i) Frictional pressure drop along the channel piPp 

(ii) entrance and exit permanent pressure losses, A Pg. 
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(i) APp occxirs aue to the shear stress on the flowing fluid 
by the walls of tne flow channel. This is given by the relation 

P = t K As, ^ ^ ( Eqn 4.7 ) 

De 2g 

where/ f''=friGtion factor 

41= segment axial length 

"•S' = water density at Tg 

y g= average channel coolant velocity 
D^= E^i^uivalent channel diamter/ and 

/ 2 

g = gravitational acclercition constant/ 9.81 m/s 

„ .p. _ flow = ^ 'a' dg.fined by 

ere / wetted Perimeter 2^r * 

(Eqn. 4.1) . 

The value of 5 depends on flow and surf (?ce conditions in 
channel. The flow conditions are characterized by Reynold's 
number , Re . 

Re = P ^e ( Eqn ,4.8) 

1.1 

where/ }j,= Coefficient of viseosity of water. 

Re should be greater than 2000 for turbulent flow that is 
desirable duexto the higher heat transport rates. In our caso/ 
237842.1. ( see p. 28 of computer printouts ). 

The second condition/ the surface condition is ch^acterized ; 

by ratio of sxxrface roughness (e ) height/ to the equivalent diameter 
. Tne flow channel exhibits smooth tube behaviour if e<t / 
the laminar sublayer thickness. Hence the cladding outer surface 
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snioothness should be so controlled that this criterion is met. 


Ihe smooth circular tube friction factor ('F . ) 

cire 

is got by using tnis coorelation widely accepted in reactor 
work; 

fcire= 0.184 * r“0’20 ^ 

The Deissler & Taylor correlation for the correction 

of circular tube data to the rod-bundle geometry has been used:- 

f 2 ^ 

vt — = -3.475 + 8.053 (^) - 4.705 (-— — ) + 0.9l62(f_)r 
^circ e 3De 

.... (Eqn. 4.10 ) 

This value of f ’ has been used in Equation 4.7. 


(ii)<lipg losses are the result of increased viscous energy 

dissipation resulting from increased turbulent activity 
at inlet and outlet. 


Pg = K 




2g 


where K = resistance coefficent for siadden contraction 
from large to small flow diameter; K=0.5 at core inlet, 

K=l,0 at outlet. 


,4.3.5 HEAT TR/J'^SFER COEFFICIENT (h) 

A two-step procedure has been obtained to calculate 'h ' 

(i) 'h ' for flow in circular tubes in found first, 

(ii) A correlation correction is applied then to obtain 'h ' 

for actual reactor coolant flow geometry. No ' -inai^'ical 

methods for predicting 'h ' for turbulent flow . . , 
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exist. However, by use of an analogy between heaib & momentum 
transport, dependence of *h' on Re Ccin be predicted. Then, 
by applying the results of a heat transfer analysis for laminar 
flow over a flat plate, an estimate of the prandtl number, P^, 
dependence on heat transfer is obtained. In this analysis, the 
Dittus-Boelter correlation has been used s- 

N^cir = 0.023 * * Pr^*^ (Eqn.4.11) ( 


where N,, . = Nusselt number for circular tube 

^circ 

Re= Renault's number (given by Eqn.4.8) 

n 0 i-i 
K 


Now, Deissler & Taylor correlation for the actual coolant 
cnannel is done, like in Eqn.4.10 


Nu 

Nu . 
cure; 


= - 3.475 + 8.053 (|-) - 4.705 (|- ) + 0.9162 ) 




^ h D 

since Nu = e 


(Eqn. 4.12) 


h = 


K 

Nu *K 
D 


(Eqn. 4.13 ) 


4.3.6 TEMPERATURE INCREASE BY THE FLOWING COOLANT 


Coolant temperature varies from a minimum at the 
inlet to a maximum at the core exit. The temperature in our 
analysis is assumed to very discretely from one segment to 
another. The bulk temperature T at any point (required in 


\ 
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Section {4,3.3 ) is determined from an energy balance between 
tbe total energy added by heat transfer/ and the energy rise 
of the coolant; i.e. 

& . ^ ^ i'-l ■ 

where, A = mass flow rate of coolant 
^1= length of the segment 

At= increase .in Bulk temperature from one segment 
to another, {Tg 2 ~ 

Thus, T-- = ^ (Eqn 4.14) 

"JL ^ “ 

in p 

for the first segment inlet temperatxire and thus is 

evaluated. For successive segments Tg 2 becomes and the 

calculation proceeds. 


4.3.7 SPECIAL CONSIDERATIONS 

Two important checks are required to ensiore proper j 
coolant conditions;- i 


(l) No Bulk-Boiling occurs anywhere in the core; 


To ensure this T^ <;^ T 


'SAT 


/ the coolant saturation 


temperature at the local pressxire. The critical location is, ^ 
of course the core exit for the central rod where bulk temperature 
is maximum and pressure the minimum. 


S3 


“Hie CHFR (CHF ratio) defined as ratio of CHF and q", should 
hence be grdater than 1 to satisfy .condition (2) above. 

4.3.8 OTHER EMPXRIC^X . GORREIATIONS 

The thermal conductivities/ viscosity/ density 
and specific heat at constant pressure are all functions .of 
temperature. As temperature in the coolant channel continue 
varies/ correlations for them would be required . 

They are (in c.g.s. system) 

Kfu(T4jr) = 0.10308 - 0.14186 E-3 * T + 0.73263E-7 * T^z 

K (T) = 0.14194 + 0.6924E-4 * T/ 

* 

for the coolant 

K (T) = 0.27656E-2 + 0.29731E-4 * T - 0.10152 E-6 * T^. 

(5?)= 0.329E-2 - 0.12479E-4 *T + 0.15296E-7 * T^; 

"S (T) = 0.7,38 + 0.6445 E-3 *T - 0.45827E-5 * T^; 

C (T) = 17,677 - 0.1104 * T, + 0.2338 E-3 * T^ ; 

p ■ 

certain other phenomena depend on local pressures. They 
are correlated as under ( p in p.s.i. ) 

Tg^^ (P) = 0.23357E3 + 0.6129E-1 * P - 0.5666E-5 # P^/ 

C(P) = 0.11312E1 - 0.78E-3 * P + 0.3312E-6 *P^ - 

0,56355 E-10 * P^z 

M(P) = 0.48141E-1 + 0.22471 E-3 * P + 0.8302E-9 *P^/ 

(C & M as in eqn 4.15 ) 
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Finally/ a series expansion has been used to evaluate Bessel’s 

functions J (x) & J- (x) and modified Bessel’s functions I (x) & 
O 1 o 

I^(x) {from ref (40)). 

4.4 RESULTS AMD C0MI4EMTS 

Analysis of core performance merchant ships 
SAVANNAH/ MUTSU and LENIN was done using the Thermal Hydraulic 
CODE Sc results obtained were within 1°C accuracy. A pellet 
diameter and clad thickness of respectively 0,858 cm & 0^.66 iran 
was chosen after an optimization study using ref (14), Results 
for various combinations of pellet dia. and clad thickness for 
a constant volume ratio of fuel & water (for satisfying Reactor 
Physics considerations) are given below; the design chosen 
(case 1) can hence be justified as being optimum:- 


Parameter 

Case 1 

Case 2 

Case 3 

Case 4 

Pellet size 

0.858 

1.000 

1.000 

0.700 

Clad thickness 

0.066 

0.066 

0.100 

0.066 

Pitch 

1.30 

1.766 

1,810 

1.262 



— 

182.41 

186.99 

130,34 

Core Power Der 

r 

*63*51 

34.42 

32.76 

67.41 

Max fuel Centr 
line Tempr, 

— 

e- 

1696.8 

1678.6 

1688.4 

1725.3 

Min. CHFR 

obtained 

2.701 

I 

1 3.334 

3.530 

2.313 

Max. biornout 

1.615 

2.946 

3,000 

1.809 

Bxxrnup^^^ 

1469 8 

16982 

16530 

10685 







































60 


The thermal Hydraulic CODE can be used for various flux 
functions Sin( z ) ^ z ,Sin( z ) & ( it -z ).Sin(ic -z ) by 

feeding an input 'flux code' 000/111 or 222 respectively. 
Some results obtained by this exercise ensue as follows 


Parameter 

Sin ( 2 ) 

z . Sin ( 2 ) 


Minimum CHFR 

3.248 

2.701 

2.797 

N 

Total F^ 

3.705 

4.367 

4.367 

Max. fuel C-Line 
tempr . ( ° C ) 

1541.4 

1696.8 

1696.8 

Max. dadoing tempr, (°C) 

335.86 

335.83 

335.88 

Flxix beyond wnich bulk 

boiling/burnout occiurs 
(times of aver?ge flux) 

1.615 

h' ' ■ ■ ■ ■" '■ ■" ‘ 

1.586 

1,586 

Coolant exit tempr. 
from hottest cnannel(°C) 

316.33 

316.92 

316.91 


Results obtained from the CODE are given in 


pages 28 to 36 of computer print-outs; flux distribution 
illustrated is Sin( z). For comparitive purposes, fuel-rod 
conditions for the hottest channel for Sin(z) & ( it -z) . Sin( it 

distribution are tabulated in pages 46 to 48 of printouts. 


’k'k'k'k'k'kikicic'kic'kick’kic 
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5.1 INTRODUCTION 


A nuclear reactor is an assembly of a proper 
size of fissionable material and other associated materials 
essential to sustain a chain reaction. Tjrie power generated is 
directly dependent on the neutron flux. It is of prominence to 
understand the behaviour of these neutrons in the reactor. 


In this chapter, we consider those reactor-core 
calculations that provide effective multiplication factor- 

and the stationary neutron-flux distribution at Steady-state 
operation. is required to establish nuclear composition and 

configuration, which satisfies criticality and control requirements. 
The steady-state flux distribution must be known to calculate 
reaction rates and power distributions that are needed for the 
thermal & shielding design of the reactor. 

Evaluation of necessitates the determination 

of the 4 factors ‘n»s»p, f and non-leakage probabilities and 
Lp . For sketching the flux distributions, 2-group 2-region theory 
has been used with water as the reflector, Heterogeneous effects 
have been considered at all stages of calculation. 

5.1.1 Neutron Cycle : 

A typical history of 100 fast fission neutrons 


and how they sustain a chain reaction is illustrated in Fig (5.1) 

All important neutron interactions are indicated in it; the 
parameters used are applicable to the present design. 



! 

i 

t 


100 

FISSION 

NEUTRONS 



'\^V'=2.5#/fission 



\ 



FIG. 5-1 CORE NEUTRON CYCLE 


I 
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A computer COPE has been generated for determining 

k and the flux distributions. The amount of 'poison' required 
err 

in control rods has also been estimated. All parameters used 
have been corrected for the operating temperature and pressure. 
Ihe inputs like core dimensions and temperatiures of fuel & 
moderator have been taken from Chapter 4 • An enrichment of, , 

m 

5,0% has been chosen (chapter 3 ). The scheme of fuel loading , 

in reactor is as per chapter 3 . 


“Ihe scheme of computations is as follows. Input 
data are enrichment, fuel-rod & reactor dimens^ns, _ref Iectqr^_^ 
thickness & average fuel & moderator temperatures. Ihe only 
variable parameter fed in is the volume fraction (of moderator) 


of Pois on, ^frac changed by interpolation till the 

critical determinant = Oj this will denote a steady-state 


reactor operation with the nett control poison of 

flux in both core and reflector are sketched for tbe above.,, 


5.2 CRITICALITY OF THE RE/iCTOR 
5.2,1 Evaluation of 

Considering the reactor core as homogeneous will 
be err«r\, ecus as neutrons may undergo several collisions within 
a fuel rod before encountering the moderator. Calculation of 
kj^£=-j|; e*pf still holdsf but differences with respect to . , , 



65 


us considerations should be noted. 


ETA ( Tl|- ) 


This gives the ratio of neutrons produced to 


that absorbed in the fuel. rod. .Since absorption cross-section 


of Oxygen is essentially zero, 

•yi ^25^f25 
^ ^a25'^ ^a28 


(Eqn 1) 


where, ' 1)25 “ 2.47 and ^ 25 ' 2^25 dependent on fuel enrichment 
(hereafter), and cross-sections corrected for average fuel 
temperature .obtained from chapter 4. 


THERMAL UTILIZATION (f: 


This gives the probability that a thermal neutron 
is ^sorbed in fbe fuel itself with respect ,to the who.le core 
i.e. ratio of neutrons absorbed in fuel only, to that absorbed 
din fuel, moderator & other structural materials. 


f - 




where, all fluxes and cross-sections are average v:,alues. Since 


0 = Cf 

^Tss ^Tm 


f = 


I 4r » V _ 
*afu fu 


Z 4 : V _ + <2 ■ V +'2' - V )t_ 

^afu fu am m ^ss ss »T 


, ..EqnC2; 


whete^-. = ^TH# is called thermal Disadvantage factor (TDF) 


(ranges, from 1.00 to 1.05). 
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This factor is peculiar of heterogeneous reactors# taking 

into consideration reduction of thermal flux in *fu^_l- rod- 
due to self-absorption. As ^^is unknown# 

f is evolved through lattice functions for cylinderical fuel 
rods:- 


1 . + siy.z) 

f ^'af u' r'u 

.... Eqn ( 3 ) 

where Ft’x) and E(y# 2 ) evolved from diffusion theory# are 
lattice functions defined by :- 

X * I (x) 

F (x) = 


E (y# 2 }) 


2 I^(x) 
2 2 

_ z -y- • 

“ 2y 


I .(z) 


O 


+ (y) 


I^(z) 


1^(9) K^^y) - k^Tz) l^(y) 


where, x= 


fu 


y = 


m 


m 


a = fuel rod radius . 

b = radius of equivalent fuel cell = 
thermal diffusion lengths, ' 


Pitcdx 

y K 


The above equations, however, are valid for b and 

inaccurate for closely packed lattices. For x, y, 
individually being ^ 0.75, the following series expansions 


accurate :- 
F(x) - ■■ ' — - 


21 


e= 1 -I? 


r (^) 


r 2 *y^'2 


Eqn (4) 


E<y,z) = 1 + i (^> ‘I ’ 


2(i-2) 


2 2 
z -y 


i»2 


. . .Eqn. (5) 


. i 
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‘f* is the most important factor that controls Control 

rod positions greatly alter it. The above treatment does not 
include control poison. In Eqn.(3)/ the value of has to be 
modified to include the effect of Boron Poison. This is done 
by considering the poison as being equally dispersed in 
moderator and a term ( gets added/ where 

macroscopic absorption (thermal) of poison & is 

the volume fraction of moderator the poison occupies. 
varies Tine^^ly with control-rod position. 

RESOM^CE .ESCAPE PROBABILITY ( p ) 

This is the probability that a fission 

28 

neutron will escape capture in resonances in U as it slows 


down to thermal energies. It is given by 


p = exp 

-'m '%m m m 


) 


Eqn. (6) 


where 0 ^ = fast f l\axes 


Here we are neglecting fast absorptions in cladding material 

a 

and Poison which are negligibly small, ^m / above# called 

the disadvantage factor for resonance neutrons is xinknown. 

So p can't be evaluated through Eqn. (6). In a similar manner 
to f/ p is evaluated through lattice functions. 

Defining 'f * , resonance utilization factor 

f = Fash absorptions in fuel resonance . ; . 

XT 

Total fast absorptions \ 
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’ I V _ - 
fu • ef f fu fu 


fu ^fu ^fu '^m'^ra 


(or) 1/fj- .= 1 + 


■fS ' V 

< Qm j 


sm , .ffl 


0 


N - I • V 0 
fu eff fu ^fu 


and from Eqn, (7) p = exp (-t; ) 

■i - 1 


Eqn. (7) 

Eqn (8) 


As ;( 




m 


w — ^ — ) is unknown, f can't be evaluated. We hence use 
^fu ^ 


the., following, lattice -functions (for cylinder ical rods ) :- 
1 


f- N 

r 


f • V 

i sm m 


fu ^eff ^fu 


* F (x) + E(y;2) ..... Eqn. (9) 


where, F(x) and E(y,z) have same meaning as Eqns. (4) & (5) 
with x,y ,2 being functions of fast diffusion lengths here.^ 
Knowing f , p is calculated using Eqn (8). In eqn, (10) i 

is taken directly from tabulated experimental values - for 
water 1.46 cm ^ ( P.232 Ref, ( 23 ) ) und the effecti^^ 

resonance integral, is well represented for cylinderical fue^ 
rods as 


Eqn. 

where A,C are constants (Ref. ibid) 
a = fuel rod radius 
"8 = fuel density. 


eff 


= A + 


. „ , nutw*! m 


( 10 ) 


(lO) 4 evaluates itself to 21, 254 barns* 






\ 


\ 
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t 28 . 


Note ; Though listed values are for pure U they are valid 
for slight fuel enrichments (upto 5%) as it does not 
significantly after the extent of resonance absorption 


in U 


28 


In Eqn. (9) the effect of Poison onjS should be 
• ^ sm 

considered*. The. effective values .of OT6ff^ .should 


be taken. 


1 


( 1 2 : ) + 

sm m 


("^2sp) 


ef f 


+ T 

sm ~sp 




Eqn (11) 


and ~ ^ 2 ) 

smeff sm sp 


FAST FISSION FACTOR (^ • ) 

All neutrons don’t succeed in escaping from fuel 
and subsequently slowing down in moderator. Some are absorbed 
in Resonances above fission threshold as well as upto the 
lower threshold of inelastic scattering in fuel# giving 
rise to fast, fissions# a favourable phenomenon, e is the ratio 
of number of neutrons having first collision in moderator . 
to the total produced in thermal fission. For determination, 
of § # the method of ^pinard# Fleishman & Soodak (SFS ) 
modified for closely-packed lattices has been 'used. The 
neutrons are classified into 3 energy-groups in fue.l with 
corresponding neutron escape probabilities without having a 
collision in fuel as P^#P 2 ^ W.^ner’s approximate 
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3 ? + 1 

it 


for i*l to 3 Eqn (12) 


where 




of fuel rod 


Surface area . 

= 2 '• for cylinder leal- rods. 

= Total Group i cross-section. 


The 2 fractions of fission, neutrons in the 2 energy groups^ 
with the first above and the second below the fission threshold 
of 1 MeV^ are given as = 0,561 and "^^2 0.439, The 

group contains neutrons scattered from the first group intc 
energies below fission threshold given by fractions:- 




V 1 ° 1£ °li. 


it 


22 


9 ‘ 0 

^ 2t 


• P = 

' 3 


Eqn (13) 


33 


3t 


NoW/ c is defined as the sum of 3 probability terms as :• 


e = 


1^1 


1 







1-(1-P2) p2 


^It^ 


y * 




1-(1-P3) P3 


Eqn (14) 
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The^, cross-section used in Eqn.lS and Eqn. 14 are taken directly 
from p. 402-406/ Ref (23). 


k. . 
inf 


Having evaluated the above is given by 


k. ^ ^ f 

inf T 


5.2.2 LEAKAGE AND 


S ince the material & geometric 


kin£-l 

buckling have to be equal for a critical system £ 


th 


taking fast leakage into consideration :- 
= ^ jnf * ^ 


eff - 1+ L 2 . b2 

th 


Eqn ( 15 ) 


In steady-state reactor operation/ ^Qff~ 1*00 

5 . 2^,3 CONTROL RE’^UIREMENT;- 

The amount of Poison required is determined 
from the CODE by trial and error method for varying values 
of such that the critical determ, inent ( refer Eqn. 63) 

is a minimum. This will be the steady— state value of Poison 
required to maintain criticality. The inbuilt reactivity 
of the core is; determined by putting ^hd setting the 

core at room temperatures. The boron Poison required is 
evaluated by choosing a slightly higher value .so that the 
reactor remains stobcritical even in the cold condition. 
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5.2.,4 TEMPERATURE AND PRESSURE EFFECTS 


(i) Ttiermal group constants 

(a) The averaged thermal absorption and fission 
cross-sections are strong functions of temperature. Taking 
the Maxwellian distribution, they vary as 


d (T) 



T ^ d a 2200 






* 


c ^ 2200 
f 


Eqn. <16) 


Eqn (17) 


where, T^ = 295,1 k T = average fuel temper atxore 

g^(T) , g^(T) are non— absorption and fission factors,^ 

respectively for deviation from Maxwellian distribution. 

For materials other than fuel, g (T) = 1.0 V T 

a 

Ihe cross-sections for all materials have been corrected as 
above using following average temper atxires 

jybderator (in core ) ........ From Chapter-4 CODE, 

Water ( in reflector) 280°^ 

Fuel temperature From Chapter-4 CODE 

Other materials in core ...... Average Coolant temper atiire 

(b) The diffusion length for water is a strong function, of 
temper atxare and pressiire. An experimental value for 
given by Deutsch ( p.92. Ref (22)) has been used; 
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2.719 

th ^ 


* ( 


T+273.1 


295.1 


0.5 


1.0614 


1+0,0614 (X,±^|Z^)0.b) 


iO.5 

) 

Eqn (19) 


where#, pr *» w.uter pressure# psi 
T = water temperat-ure# °C 

'B (<];') = water density at p^ and T 


(c) For water# 


2^ o<; N# nximber density# which in turn is 
<K "§(T) J '2:, ' 


As# 'S(T) = 0.93796 + 0.64448E-3*T - 0.458271E-3 *T^ gm/cc 

Eqn, (20) 

where# T= water temper atxire (for 2000psi pressurised 


water ) 


So# 2 (T) ='5(T) * (^)0.5 ^2 (Tq) 


^ 

O 

Next# to evaluate diffusion constant ! D for water# we use 

^ ^ 

the_ familiar equation 


D = *y 

fli -^a 


^ Eqn (21) 

We see that the D calculated in the above manner is 


reasonably accurate. 


\ 
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(i±) Fast group constants 


(a) The neutron age# i , ixi water is not 
that much affected by temperature as pressure. ic ; 

As temperature & pressure increases# density decreases 
thus increasing neutron age in water. We arrive at : 


% (T,pr ) = 


0 


Eqn. (22) 




where 


2 t 

= 27 cm X at atmosiiheric pressxire 

= 1 gm/cc j and temperature of 22.1°C 


0 


(b) Ihe fast neutron age in core is dependent on temperature 
in an indirect manner. A temperature increase results in a 
decrease of thus increasing the age. This has been elaborated 
in the subsection that follows. 

t* «■ 


(iii) Heterogeneous Core constants: 

(a ) In a heterogeneous lattice# the 
average absorption cross-section '2. * is 

3.0 

as afu fu ap^ ^ ' ass ss 

.. Eqn (23) 


and 


■^ = 2 - V + ^ 

t^r trfu fu -ferm 


^"trss ^ss 


- Eqn. (24) 


*trp being relatively very small may be neglected. 


* Hereafter subscript 's'. .stands for thermal/slow neutron 

parameters. Subscript 'f ' stands 'for, ^Fast neutron parameters. 



(b) The extrapolation distance (E.D. ) of the core is then 


E.D. 


0.71 




^'tr 


Eqn(25) 


(c) Since average diffusion coefficient, D , r.f 'Loctice varies 

o 

inversely as 2^^ / it follows from Eqn.(24) that 



V 

D 


fu 




V 


in 



+ 


V 

SS 


D 

SS 


(Dt) 

s 


1 


fu 

, 

V 

. SS 

fu 


H — 

SS 


Eqn ( 26 ) 


(a) Now hs t the thermal diffusion length of core, is determined 


from Eqn (26) 5<^(23) 
D. 




2^s 


Eqn (27) 


(e)^‘-The neutron age in the core is determined from experimental 
values for age of a |^H20 mixture (p.206 Ref . (23) ). In the 
above case, a change in temperature ( & pressxore ) has the 
effect of reducing the atom density of water relative to 
fuel which is equivalent to increasing the metal- to- water 
ratio in voltome. 

As “S is given by Eqn. 20, the ratio of U-H 2 O mixture, say, 

X = Vol.of U * 

Vol. of H 2 O ^(T) 

(23) 

Using Fig 6.4 of ref , we arrive at the following correlation 
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t; 19.8 + 12.5 * X - 5.0 * X? ^ for X -^0.9 
= 27.0 + 6.27 * (X_0.9 ) ^ for X> 0.9 


(f) In Eqn (26) D is taken from Eqn (21) but D . D and D are 
^ m ss* p fu 

still unknown. These are individually calculated as 


1 

“ r (:% 1^)) 


Eqn (28) 


(g) Finally / the fast diffusion constant/ in core since 
hydrogen scattering cross-section varies significantly with neutron 
energy/ an artifice is used for core transport cross-section. The:v.. 

using j 

V Eqn (29) 

"-af = ' 

m 

and Eqn. (28)/ is calculated. All the above calculations are 
shown in Apx^endix A5-1. 


5.2.5 HETEROGENEOUS VERSUS HOMOGENEOUS 

We are now in a position to comprehend what 
effects would weuhave incurred had a homogeneous treatment of 
the reactor been done. As '1'^ "/“ij f (hetero) ^ f (homo) due to 
the self-shielding effect of the fuel rods. However/ p (hetero)^ 
p(homo ) as the above effect reduces the number of neutrons 
that are parasitically absorbed by fuel resonances. Thus/ decrease 
in f is more than offset by increase in p . So/ (f.p) hetero^ 

(fp) nomo. 
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The fission neutrons pass through a region of pure fuel where they 
are likely to induce fast fission, before emanating to the 
moderator. Hence, ( e) hetero ^ ( s ) homo. This effect is not 
that pronounced as the other factors above. ^ does not, however, 
change. 

In view of the above inequalities, for 

heterogeneous lattices is greater than its homogeneous equivajjit, 

A homogeneous treatment, though much simpler would consequently 
give erroneous results. 

3. FLUX DISTRIBUTION 

The thermal reactor consists of a heterogeneous 
multiplying core surrounded by a non-multiplying water reflector. 
Since the core consists of a closely-packed square lattice of 
fuel rods, it can be considered homogeneous for determining flux 
distributions in core Sc reflector. Despite its limitations, 2-group 
apx-roximatioiQ has been used, as it is the simplest that will give 
reasonably accirrate results for multi-region calculations. 

The diffusion equations for fast & slow neutrons 

JtS" given by :- 

- Ij-f’^f + ^Inf =0 Eqn (30) 

P-2’rfh = 


0 


. . . . Eqn ( 31) 
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where, all symbols have their usual meanings , 


= thermal absorption cross- sect ion 
= slowing-down (fast removal) cross-section. 

In core, 1 in reflector 0.0 

The solution of Eqns (30), (31) are solutions of wave equations 

7 ^ = 0 

7 ^ = 0 

This implies the coupling of the solutions by the condition the. 
coupling- coefficient, S , between fast & slow fluxes is 


S = 


+ D B^ 




~’as 


"V 


Eqn(32) 


rf 


The criticality condition should be satisfied; 
k 


i.e. k = 

elf 


inf 


Eqn. (33) 


( 1 + 

where x = 1/L 
s s 

We shall use subscripts '1' for core and '2* for reflector. 
Eqn (33) defines 2 values for fig. 5. 2) 

,2 


2 • ' vf --^ ■ ^ 2 


B‘ 


- + 




s 


Eqn. (3&) 
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B. 


( xi^„) - V^( + x^_)^ + 4 Xt „ x,^ . 


Uf ""Is' 


^if-: is; 


.Ic . 1£ inf 
Eqn. (35; 


2 2s 


n' 2 _ 2 

®2 “ ^ 2f 


.... Eqn . ( 36 ) 

Eqn. (37) 


Considering wave-equations are separable in botii core & reflector/ 


0 


Cos (- 


% 


H + 2RS 


-) * F(r) Eqn. (38) 

D, 


R^-R 

where/ RS = reflector saving* = (■ ’ - ^ — )- L^k ( 1-exp (-rp — ^ ) 

O- - •z^s 


2s 


2 s 


F(r)- is the solution of the equation 


d F 




dr 


since B has 2 values in both core & reflector (given by Eqns(34) 
to (37) ) / to define F(r) conpletely/ a total of 4 equations will 
have to be satisfied s- 

+ 1 - d F^^ ^2 = 0 Eqn. (39) 


dr 


-.2„ 

d F 


If 


dr 


dr 


+ ^ ^If - F ^ = 0 Eqn (40) 

r dr 


12 
d F 


dr 


. d F 

. ^ 2f 

2 r dr 


.2f 


1^) 


F2f = 0 


. Eqn (41) 


* From p.69, ref (24) 
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^ ^2f -I- 1 ^ • F =0 Eqn (42) 

.2 r dr s 

dr 


Identical 4 more equations are there for 

m, and are all +ve constants 

2 2 

..2 ^ 

1 ^H+ 2RS ^ * 

2 2 

2 _ /_ 2 L __4 

" H +2PS‘^ 

9 9 71 ^ 

“ ®2 ^ H + 2RS^ 

2 ^ 

“ ®2 H +2RS 


and F^ as well. 
2s 

Eqn (43) 

Eqn (44) 

Eqn (45) 

Eqn. (46) 


The solutions of Eqns (39) to (42) are Bessel's functions. The 
complete solutions are as under 
In core 

^Is^^^ = A-’J (Ir) + C'l (mr) Eqn. (47) 

0 0 

F^^(r) = A-S^-J^dr) + C€^» Iq (mr) Eqn (48) 

In reflector 

F^g(r) = E*Iq( M’^r ) + G-Kq ( ^^ ) + H'lQ(l^£r) +M.]<Q(n^r) 

Eqn. (49) 

F2£(r) = S^\( H-Iq ( ^^r ) + M-kQ(P-£r)) 


Eqn. (50) 
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wiaere s given by Eqn. (32 ) are 

,2 


= 




X, + D • B' 
las Is 1 


P-^ 


Eqn (51) 
. Eqn (52) 


Irf 


S 

2 



= 0.0 


V ^ 

•^'las 4- 25^^ 

^2rf 



Eqn (53) 
Eqn (54) 


Eqns. (47) to (50) contain 6 unknown constants A, C, E, G, H 
and M. These can be determined by the following boundary 
conditions 




^2s = 0 


Flf(Rl) = 




If 


dr 


= D 


2f 


dr 




dr 




dr 


Eqn (55) 

Eqn (56) 

Eqn (57) 

Eqn (58) 

.... Eqn (59) 


Eqn. (60) 
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Applying Eqns (55), (56) to (49), (50) gives 


E 


^0 ^ 


0 


Eqn. (61) 


H 


M 


V~ 


V2) 


Eqn. (62) 


using Eqns (57) to (62) in (47) to (50), we arrive at the 


following 


A 

C 

Yll 

Y12 

Y21 

Y22 

Y31 

Y32 

Y41 

Y42 


G M 


Y13 

Y14 

■ 

r— — , 

0 

Y23 

Y24 


0 

Y33 

Y34 


1 0 . 

Y43 

Y44 


0 


— 


~ -1 


.Eqn. (63) 


wnere, Y. . 's above are elements of the critical determinant 
•IJ 

as given in the CODE ( line nos 23900 to 25900 of main 
program ), /' .■ ut;,. 1 1 ..i:i ..j 5 'x X‘.i.l--c v*. j..: tii*^ '-‘i’ ■■■'' 

C’ Yhe pciti'-ii deterrnir, '-nt 

; only then will A,C,G & M will have non- 
trivial values as the 4 equations are homogeneous. 

Next^ the first 3 equations of Eqn (63) are taken and constants 
G, M and G evaluated with respect to A p,sing Crammer's rulei, 


j. . 


given by 


Eqn (63) = 0 
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Finally/ A is determined by the following method. Since the 
average thermal flax in the core is given by the relation 


+H/2 


TCR^H 


< A * Jq ( Ir ) + C ' Jq ( mr ) ^ • Co s ^ 


-H/2 0 


2 n rdrdz 


from A 3.1- 3 ^ this works out 





I^(mR^) 

mR^ 



) I 

(H4-2RS)*2 / I 


Hence/ 

A. ^ * S3 — 

?sav + 2( 2<H-2RS)Sin^5^j } 

IR^ ^ mRj 

Eqn (64) 


It should be recalled that { ■^ ) in Eqn, (64) is known and we 
obtain 'A’ as the number of times of From this/ the other 

5 constants can be evaluated and since all constants in 
Ec^s (47 ) to (50) are now known the radial fl\ix fli notions in 
core & reflector can be easily sketched and is shown in Fig. 5,3 




5.4 RESULTS 


A simple one group calculation of versus nxamber 


.28 


den. ratio of U & water as per p.306/ ref (23)/ gives a maximtom 
at a ratio of 0.156 for the nuclear parameters chosen; this 
implies a pitch of 2.12 cm for a fuel rod diameter of 0.858 cm ai 

3 

a core power density of 24.52 MW/m . This is clearly undesirable 
for the compactness desired. However/ a pitch of 1,30cm chosen/ 
though resulting in lesser moderation is a compromise/ with the 
operating point falling only 8% from the peak. 

Core criticality has been achieved at k^^£=0, 992945 
with the aid of a 10.99 cm thick water reflector. The various 
intermediate parameters/ criticality determinant & others are she 
in file CHECK-OUT ( p.39 of computer print-outs). Other nuclear 
parameters like cross-sections / etc. are displayed in file K.OUT 
(p, 37-38 of Print-outs). Fig 5.3 maps the radial flux variation; 
a rise in thermal flux at core edge due to moderation by the wate 
reflector is clearly seen. ! 

Results of reactor control calculations (refer chap 7 
evolves a volumetric Poison fraction of 0.92148E-3 for a clean ; 
reactor with a maximum of 2.35E-3 (with all rods in) to drive 
the reactor subcritical by 5% /^k/k in the cold condition. File ‘ 
PSN.OUT (p.40 of print-outs) lists some parameters discussed ; 
later in Chapter 7. ! 


’:kic:kicie'k'kic'k-k'kii:ic 
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REACTOR SHIEUDING 

RESUME 

6 . 1 INTRODUCTION 

6.1.1 Types of Shielding 

6.1.2 Shielding Requirement a 

6.2 SALIENT JFSATURES 

6.2^1 aiielding Design Characteristics 

6.2.2 Dose Rates Permissible 

6.2.3 Optimization Techniques 

6.3 RADIATION ATTENUATION IN SHIELD 

6,3. 1, Nautron attenuation calculations 

(i) Fast neutrons 

(ii) Thermal neutrons 

6.3.2 Gamma Ray Attenuation calculations 

(i) Core gammas 

(ii) Capttire gammas 

6.3.3 Emissions from reactor effluents 

6.3.4 Total Dose -rates 

6.3.5 WeiQ^ta of. shielding materials 


6.4 


COMMENTS 



INTRODUCTION 


Shielding of a marine 


reaetor is ef^faT'Ijgreafcer '’i%'ortance than its landbaeed 
counterparts.^ Here^ the peculiar space/ weight and dose rate 
limitations require an optimization study concerning the 
shielding materials and the differential shielding thicknesses 
on the btasis of dose rate mapping. Reactor shielding is 
essential to t- 

(a) Not only protect the personnel working in close proximity 
to the reactor system but also to 

(b) Reduce radiation levels in the vicinity of instrumentation 
and other equipment that may be damaged by radiation; 

(c) Reduce the. induced radioactivity of reactor components 
that may need to be serviced during shutdown (like primary 
loop components) 

(d) Prevent excessive radiation induced heating of tl^e 
pressure vessel (PV), 

(e) Limit the neutron erabr it t lament cf the PV, 

' The marine reactor shielding design would require 

a careful study, of propagation of radiation throu^ matter .. 
(i.e. shield material), system analysis of activity production 
and release in reactor effluents* It also requirea other 
aspects of the ship. & reactor design such as general plant 
layout ^ structiaral studies (not done here), Ihe aim ie to 
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utilize the large amounts of materials needed for shielding 
also for structural purposes to optimize on weight. For 
analysis a ^MARINE REACTOR SHIELDING COE® has been developed 
which has been referred to throughout . 

6.1.1 TYPES OF SHIELDIfJG 

There are 3 Key factors in shielding design- 
cost/ shield vol\ime & shield weight. Any one of these may 
decide, the type of shielding required. There are three types 
of shielding:- 

(i) .Unit Shield - that completely sxarrounds the reactor / like 

in land-based reactors. - 

(ii) Compartment Shield - that shields the reactor proper but 

not all its components like in 
mobile reactors, 

(iii) Shadow Shield- that shields only the crew and important 

instrumentation, like in aircraft or 
space reactors. 

'i 

6.1<2 REACTOR SHIEl^DING REC^IREMENTS 

The oontrol of radiation from a ^lip-board 
reactor presents, in particular, the 2 following problems 
of design ;- 

-shielding from, neutron & gamma sources in a compact geometry, 

layqpl^, 

-management of safe release of gaseous activity to satisfy 


\ 
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the requirements of rapid access to the reactor containment 
after shutdown. Among these, only the former has been 
studied to some depth. 

The scheme of reactor shielding revolves 
primarily aroiond these 5 considerations 

(a) Slowing down the fast neutrons 

(b) Capture of slowed down on’. 

(c) Determining the soiarce of the secondary gamma radiations 
formed as a result of various neutron interactions in 
shield materials. 

(d) Attenuation of all forms of gamma radiation. 

(e) Determining the total dose rate of radiation emerging 
from the shield materials. 

Let us examine the abov^^ points; for high 
energy neutrons, use of elements of moderate or high mass 
nximber is effective, as they are able to reduce their energies 
through inelastic scattering. Since elements of low mass 
n\imbers are the best moderators, in the ,form. of wat^ is 
used as a shield constituent for slowing down intermediate ■ 
neutrons t Hence a combination of a moderately heavy. 

i 

element like steel with water is used* 

As a general rule, most of the f-rays 
originating from the core get absorbed in the shield region 
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nearer to the core itself. It is the secondary ¥-,rays^ 
however/ ( emitted due to thermal neutron capture within 
the shield) along with a tiny fraction of the core gammas 
that control the shield thickness required. The neutron 
attenuating material to an extent also serve this piirpose. 
However, as y-ray relaxation lengths are inversely proportional 
to material density, additional material of higher, density 
(like Lead) is incorporated to attenuate the higher energy 
f- rays. 

A word about the 'thermal shield'. This is basically 
to safeguard the PV, High gamma flux from the core leads to 
very high heat generation rates in any struct.xiral elements 
placed there. If it were the PV, then the resulting large 
temperatures & thermal stresses when added to the . already 
present high mechanical (pressure) stresses, would bring the 
total stress on it (PV) to an unacceptably high level. Hence, 
a thermal shield consisting of layers of stainless steel 
and water is placed between the core & PV, The arrangement 
is as shown in Fig. 6,1. Ihe steel has a high melting point 
& excellent thermal conductivity. Energy is absorbed by 
attenuating "f-radiation and inelastic scattering, of fast 
neutrons. Ihis heat is immense & la rerroved by the, coolant 
flowing into the reactor core through these shield layers, 
thus contributing to the nett available energy, from core. 




!fe 


Fig. 6,1 


Thermal shield in the PWI^. 
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Ihe core- ^ energy is attenuated from ^ 10^"^ MeV/cm^-a 
9 2 

to 10 MeV/cm and the fast neutron fliax from 

0.28 X 10^^ #/cm^-s to 0.16 x 10^ # /cm^-s before 

impinging onto the PV, Thus high thermal stresses and neutron 
induced embrittlement ( by high-energy neutrons) is effectively 
reduced. 

6.2 SALIENT FEATURES 

6.2.1 SHIELDING DESIGN CHARACTERISTICS 

The scheme used is as follows;- 

(a) Thermal shield: stainless steel <33 347) and water 

alternately used. 

(b) Pressure Vessel: made of carbon steely as it is required 

to withstand high mechanical stresses. 

(c) Primary shielding; also- called Neutron Shield Tank(NST) 

surrounds the PV-Made. of carbon ^beel 
Sc water alternately# it has imxx>rtant 
mechanical support functions besides 
ehie Icing# like supporting the entire 
PV#c:RI» 1# Cooling ducts# part of 
primary piping# etc and also to resist 
shock load acting on the ship. The 
concentric layers of steel and watef 
• makes an optimal shield for both 
neutrons 6c)f -rays. Lead is not used 
here due to its lower melting point 
& poorer mechanical properties. 



(d) Secondary Shielding - It sxirrounds the containment vessel 

& the complete primary loop (steam 
generator, associated piping, primary 
, pump, etc.)* It attenuates to the allowed 

dose, radiations coming from coolant 
loops ( reactor effluents)- mainly 

# r 1 A 

)f-rays from activation) and those 

escaping from primary shield. It's made 
up of one layer each of Iiead Sc Polyeth-'-ne. 
Concrete due to its much larger bulk for 
the same effectiveness has not been used. 
Study reveals that use of Boron based 
alloys here do not substantially aid 
in reducing the shield weight/bulk. 

Table 6.2 shows the thicknesses of the 
above shield materials evolved from the 
shielding CODE. 

6.2.2 PERMISSIBLE DOSE RATES 

The occupational individual accummulated dose 

is limited by the well known equation D = S‘(N-18) rems, where 

N is age of the person. This implies a maximum dose rate of 

2.5 mrem/hr. For personnel working in nuclear sites the* limitation 

however is l-.S rem/year or 0.75 mrem/hr. This is the design 

consideration in our shielding calculations. However a total 
(29) 

dose of 12 rem under emergency conditions has been allowed 

with doses^excess of 5 rems being adjusted over the next 5 years. 
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6.2.3 OPTIMIZATION TECHNIQUES 

The search for the optimum shield consists in 

evaluating the weight of various possible shield designs giving 

<28) 

same dose rates at specified points. B. Chinaglia , has 
sketched weight variation W for the optimum configuration, as 
a function of neutron (D^) to total dose Ibf) i.e. 

D^/(D^+ DY). Ihe curve proy^es that the minimum weight is found 
for D^/(D^4 -dY) = 0.10. This gives the lead & polyethene 
thicknesses as mentioned. 

The optimim primary shield with respect to weight 
is composed by laminations of heavy metal located at about 
the middle of NST. No strong difference between the use of iron 
only or iron and lead is found. Hence only carbon steel is used. 

6.3 RADIATION ATTENUATION IN SHIELD 

The radiation from core includes a , p 
and iT-rays, neutrons of various energies, fission fragments, 
protons & a-particles from (n,p) & (n,oc ) reactions and 
neutrinos, cmly neutrons & gamma rays are considered for 
shielding design as they are most penetrating. The primary 
and secondary radiations from the core have been briefly | summed 
up in Fig. (6-3). We shall consider neutrons SclT-rays 
separately. 




FjQ.5.a PRIMARY & SECONDARY RADIATIONS ‘SCHEMATIC . 
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6.3,1 NEUTRON ATTENUATION CALCULATIONS 

Ilie theory of neutron penetratii^n through thick 
layers is different from ordinary theories of neutron diffusion 
used for criticality calculations where the interest is in the 
behaviour of average particles. But here, the interest is in 
those exceptional high-energy particles that do not undergo any 
interaction c throughout the snield materials. The average neutrons 
that do the normal things are not very important as they stand 
no chance of getting f^r into the shield. 

Neutrons are absorbed in the shield in 2 distinct 
stages. First, high energy neutrons are slowed down by elastic/ 
inelastic collisions; these slowed-down neutrons are then 
captured in shield materials easily as capture cross-secctions 
are high at low energies. We shall now consider these 2 stages 
separately. 

(1) FAST NEUTRONS 


Theory of attenuation 

These attenuate as per the following 3 forms 

(a) Elastic Scattering- Here, neutron K.E. is lost by collision 

with target nuclei. Dependent on scattering 

nuclei mass, it is the most important 

form of energy loss for light nucleg. like 
where a single collision 
may result in the removal of the neutron. 

From l-eV to lOKeV, this form is dominant. 
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(b) Shadow Scattering: This occurs at energies greater than 

5 to 10 MfciV wherein the neutron acts as a 
wave and is diffracted about the target 
nucleus. It occur es only when neutron wave- 
-length is lesser than nuclear radius. At 
larg'e neutron energies/ its cross section 
approached f\.R , geometrical nuclear cross- 
section. 


(c) Inelastic Scattering: Its importance is that the incident 

neutron is strongly degraded since it loses 
the recoil energy of target nucleus plus 

, r 

energy equivalence of inelastic v -rays 
emitted. This is very important for heavy 
nuclei where allowable energy levels are 
nximerous and smaller^^^istances between ground 
state €t first energy level. Its threshold 
varies from 6 MeV for -to few tenths- 
of-gji' MeV for Uranium. At high neutron 
energies, its cross-section also approaches 
riR^. 


In marine reactor design, we are interested in 
minimum weight. Limiting values of both (b) and (c) above approach 
. Hence total cross-section (*?) ) approaches 2^^ at high , 
neutron energies. As nuclear radius R=1.5 x lO”^^ * where 

A is mass number of the nucleus, the total mass cross-section > 
becomes : 

6.022x10^-^ (2«R ) 

= A 


0.0852 

.1/3 


2 / 

cm /g 


1 



which means light elements like give minimum shielding weights^ 


justifying what was stated previously, 
Ecaiations used 


Neutron fl\ix follows the exponential attenuation law 

0^iy:) = 0^' . exp (-x.-.x) 
where x = distance from source 
0^ = soxjrce strength 

^ = effective removal macroscopic cross-section 
For 2 to 12 MeV neutrons 

0 = 0,35 *1 barns 

of material with mass no. A. This is used to evaluated of 
unlisted materials like steel (refer Appendix A6.1) 

Applying inverse square law for spherical attenuation 
and considering 0 £(R) as flux at core edge, 

0^'(x) = 0.^(R), exp (-Xr ).R^/r^ for r^ R 
As more than one layer is present, the exponential term for 
'n' layers gets modified to 




where 


^ 1 = cross section , 

thickness, of layer i 


The attenuation prior to the P.V, wall occurs 
in the water and steel of thermal shield region. Experimental 
Data (p. 153/(26) ) for this region show that when build-up is 
considered, 
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point kernel can j-e represented tv the sum of 2 exponential 
terms. For a large planar source of core/ flux distribution 
in X cms of water/ (x) , is 

0^ (x) = 0o C A*E^(ax) + (l-A)*E^(cx ) 

where/ 0o = inciaont fl\rx; 

E^( ) = incomplete error functiv-ns 

(refer App. A6.2 ) 

A = 0.89 

a = 0.129 / g= 0.091 (constants.) 

The attenuation of fast neutrons thro ugh z cms. of steel 

I 

followed by x cms of water is then well approximated by 


0^(z,x) = 0^(x) * ) 


Where 2= removal cross-sec for steel. Flux thus 




obtained is multiplied by 
stherical geometry. 


R" 


to account for the 


(R+z+x) ' 


Method used 

The method used in SHIELDING CODE is 'Diffusion 
Removal Method'. Here/ the basic assumptions are 

(1) The penetrating component of source neutrons consists of 
high energy neutrons that suffer only a small energy loss 
and uncollided neutrons. 

/ 

(2) Neutrons suffering large energy loss are regarued as being 
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removed from the beam. 

(3) The removed neutrons are degraded in energy according to 
diffusion theory and don't travel significantly from point 
of removal. 


The equations explained above have been used in 
the code to obtain the fast flux at every layer. The fliix 
attenuation through the shield material is sketched in Fig. (6.4). 

(ii) thermal neutrons 


The source of thermal neutrons are 


(a) thermal flux emanating from the core. 

(b) neutrons removed from the ^ j--, 

fast-flux beam. These removed 

neutrons / forming local soxxrces can be adequately 
described by diffusion theory. This source intensity ig 
* 21 .In any layer a distance L from core 

centre , 


= 0^ (L) * 1? /r^ * exp (-2(r-L)} 

where = fast flux at inner face of the layer 

tand r ^ L. 

These neutrons can then be introduced into the highest energy 
group of an appropriate set of multigroup diffusion equations 
to calculate the lower energy neutron distributions. Here# * 


only 2 group . . . 
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theory has been assxmed. Hence# all above neutrons are assumed 
thermalized, (This only results in slightly higher dose rates). 

If 0 (r) is the slow flux at any radius' r# 
s 


D'V^ 0 (r) - T *. 0^ir) 0^(r)»2 = 0 


/ ' ■ i 

le^age term absorption .local source term, 
term 


where D = diffusion constant of layer 
s 

2: “ thermal absorption cross-section, 

as ^ 


(or) .0^ ~ 0^ir) + Q 


0 


where 


# = V- 


H /D 
( as^ s 


Q = 0p(L)-.L^ 


D 


s 


Evaluating for = ”2'"’^ (r^d ) (spherical coordinates# 

^ ^ considering only radial 

variation) 


0^(r) 


Ae 


-cr 


Qe 


-Eh 


2 cr 2 cr 


(c.-S). Bj<y)^cr 


( refer Appendix A6-3) 
where ^ x = (c-'£)r ^ y= x/ (l+2c/{c-2»’)) 


As 


«cr 


cr 


<< e- we modify as 
-cr 


0s (r) =— • 


Ae 


2cr 


+ O e ^ (c- jE ) . E- (y) 

2cr ^ 


. . . eqn ( 1 ) 



This equation is difficult to evaluate. As the effect of Inverse 
square law is much lesser than exponential attenuation within 
a layer/ the distance attenuation is neglected within the layer 
for 0 (r) evaluation. This is a resonable approximation If layer 
thicJcness 


Now diffusion equation redqses to 

^ (r) - c^ 0 (r) + = O 

where/ Q = H e . ^ 

^s f 

The solution of this equation (refer Appendix A6.3) 
is 


Pf3(r) 


Ae""^^ ■ - r 

Qe 

r ^ “r 


( 




(c^-I? )2 


which is easily solvable; flux thus obtained is multiplied by 
2 2 

L /r for distance attenuation. It must be mentioned that this 
method gives a slightly higher value for i2f^(r) than Equation (l), 

o 

The thermal flux within every layer is thus 
obtained & has been mapped in Fig, 6.4 . These flux values will 
be required later to evaluate the capture gamma-ray sources 
(local) within the snield material. 
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6.3.2 GAMMA RAY ATTENUATION CALCULATIONS 
'HiEORY OF attenuation s 

Depending on the 1/- energy and absorber mass- 
number, the 3 absorption processes of interest are photoelectric 
effect, Compton scattering and pair production ( in MeV region). 
Secondary effects such as induced reactions, cohtront scattering 
and Raman scattering have not been considered as their effects 
are neyligible. For illustration, V" -ray cross-section is 
shown as a function of energy for Iron in Fig. 6.5. 

(a) Photoelectric effect : Most important absorption at low 

energies «1 MeV) Energy is transferred to orbital electrons 

through Einsteins equation; 

2 

Esht=w +%mv 

\ \ \ 

Total energy work JC-S-of elctron 
function 

It is a true absorption process as the V^-ray gives up all its 

energy in one interaction. The cross-section for this interaction 

is proportional to z'^/ where Z= Atomic nximber, i>=V-ray energy 

in MeV. 

& n =■ 3 to 5 depending on ray energy. 

(b) Compton Scattering : Occurs at high energies. Here, y -rays 
undergo billiard-ball type collision with electrons, 
according to compton's equation: 

E= Eo 

1+ ^ ( 1-Cos O) 

0.511 


where O = scattering angle 
Eo= incident^- energy 
E = ScatteredV«nergy 




gamma energy (MeV) 


FIG. 6'5- GAMMA 'RAY CROSS 'SECTION FOR 
IRON."" 


Fig .5.3 , ref (26) 
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It is nor a pure absorption process as fractional energy lost 
depends on incident angle. Its cross- sect ion is proportional 
to Z/E ( nearly constant for all elements (except H 2 ) where 
Z/E =0.5;) . 

(c) Pair production ; At higher V - energies, this startling 
phenomeiion occurs wherein under the effect of the coulomb field 
of the nuclues, a high y)- energy (> 1.022 MeV) can give rise 
to an electron pair-positr on & electron, with any excess energy 
appearing as particle kinetic energy. Positron always gets 
annnilated by combinatien with an electron giving 2,0.511 MeV 
V-rays. Its cross-section is proportional to Z*(E-1.022). 

It is a true absorption process as annhilation radiation is soft, 
emitted isotropically and attenuated easily. 

The above 3 process are all added together to 
give total if -macroscopic cross-section . For materials like 
stainless steel containing rrrare than one nuclear species, the 
individual are added together to give the total. 

This is snown in Appendix 6.1(b). From the table of mass- 
absorption cross-sections for ][ ' s (p-762, (21)), it is clear that 
except H^, the most efficient f - absorbers are heavy elements 
like lead ( & iron, to a lesser extent) which have high cross- 
sections for both photoelectric absorption and piiir production. 
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(i;) CORE GAI'MAS 

The reactor core exhibits a complex V" ray spectriain. 
For shielding calculations it can be summed up as follows ?- 


TASLE 6-3 


V-RAY SOURCE 

J PHOTON 

t 

ENERGY 

- (MeV) 


% 

I 1 

2 

4 

6 . 

8 

Prompt fission MeV/Fission 

i! 

i 3.45 

3.09 

1.04 

0.26 


Fission product MeV/Fission 

5.16 

_J 

1.74 

1 

0.32 



- 

Total 

j 8.61 

CO 

• 

1.36 

0.26 

! 1 

- 


Besides these/ capture V" -rays are also being continually emitted 
through thermal neutron capture in water ^ cladding material /etc. 
All this is summed up as in Apprendix A6,4. The core has been 
considered as a large uniformly distributed source of volumetric 
density (V^^) dictated by thermal hydraulic design. As 1 watt 
power is 3.1x10^^ f issions/cm^-s considering 200 MeV energy 
emitted per fission, the volumetric soiirce strength would be 

i 

* in 

= 3.1 X lO-^^- Mev/cm -s 

where, V'j^'s are chosen energy intervals from 
1 to 8 MeV . 

since the core absorbs y-rays emitted by it, this self- absorption 
phenomenon has bean considered. This reduces resultant radiation 
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level in the core outer surface by F . times / where F . is 

pi pi 

self-absorption function at the ir -energy . 

„ Radiation intensity with self- ab sorp t ion 

Sp - 

Radiation intensity without self- absorption 
For a cylindrical source viewed from the side/ 


^Pi = 


% 


R 


where/ R = core radius 

2^ =V -absorption cross-section at 
\ -energy Eih.. 


If dose rates are to be considered at distances >)>• R/ then the 
reactor can be considered as a point isotropic source of strength 

Vi 

= fay X V fk Fpi MeV/jj ^ where V = core volume/ for 

V -energy Ihis is the source strength used in the CODE, 

In the CODE/ the point kernel method has been 
employed where the reactor is considered as a source as above. 
Since V-rays follow the exponential attenuation laW/ in 
spherical coordinates/ the /-flux will be :- 


0 (r) = Sj • exp(-rr) 


where r= distance from source centre. 


4 t 

For attenuation through different layers of thickness t^^ and 
V* absorption cross sections the exponential term becomes 


- b 


all i 
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From the main beam of \’s, tne photons being scattered 
out - of any particular narrow beam are replaced to some extent 
by particles being scattered- in from .adjucent beams. This 
results in an increase in particle flux- density depending on 
if-energy & material and *'b". Hence a factor , isotropic 

point-source dose-rate build-up factor (> 1 ) is included in the 
f 1 u>;fequat ion . 

0 (r,Ey. ) = ^ exp(-b(B>l)) 

415 : r 

when the total dose is required to be calculated , the effect 
of every By. is to be summed up. This equation has been used to 
obtain the effect of various energy fluxes at different layers. 

(ii) CAPTURE GAMMAS 

These rays emitted throughout the shield, are directly 
dependent on the thermal neutron fluxes. From this, and a 
knowledge of material densities & cross-section, the secondary 
gamma source-rate may be calculated thus:- 

Source (i, r) = ^ 0s (r) photons/cm^-s 

oiS 

at layer i, a distance r from core centre. In most cases, a 
spectrum of Y-rays are emitted. These have to separately be 


considered thus. 



Ill 


Source (i, r, ^2fs(r) 

Here, again, assumed to vary in discrete steps 

from 1 to 10 MeV, 

Since V -rays are generated throughout the shield, 
consideration must be given to the varying shield thickness 
available to them. So, each layer of the snield has been divided 
into a series of concentric rings of thickness ^ 1,0 cm and each 
of these are Consiaered as separate sources, Tne total radiation, 
outside the shield is obtained by summing the radiation obtained 
from each of these rings from all layers, plus that obtained from 
the core itself, following the Point-Kernel technique (described 
unaer "Core Gammas") of \''-ray attenuation, for each energy 
interval, ^ This procedure has been adopted in the CODE, 

6.3,3 EMISSIONS FROM REACTOR EFFLUENTS 

The use of light water as coolant poses far 
lesser problem than other liquid coolant in this field. The major 
reaction of interest is « 

(n, p) activatiun cross section 2.0E-5 barns ^ 

17 

N emits gamma rays of average energy 6.2MeV-The total activation 
energy evaluates as:- 

= 8.2925 E~6 0^ ’V Vfj 20 ,eff * exp (-0.09367 T ) MeV/s 
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where / = virgin fission neutron flux 



effective water volume to be considerea for 
shielding purposes. 


r = effective time constant for transport. 


In the abscence of Vn^O 2 nd 5 and shape & cim4.nsions of all 
primary components (extra-core) required to calculate equivalent 
effective flxox for snielding/ this emission has ret been 
consiiered for reactor shielding here. However, since its magnitude 
is small, it does not introduce great errors in dose 

calculations. 


6.3.4 TOTAL DOSE RATES 

Once the neutron and gamma fluxes are known at 
all interfaces, it is easy to evaluate the total effective dose 
rates. For this, the flxxx for unit tissue dose (l mrem/hr) is 


used as simple conversion factors (CF) . Taken from P.7-66 of 
Ref ( 21 ) they are ;- 

Table, Fluxes foY- tissue d-ose vAss s 

_ ^ 

Particle j fltix for 1 mrem/hr dose. 

j 

J 2 

(a) Thermal neutrons J 1500 #/cm -s 

, . < 2 

(b) Fast neutrons I 7. 333# /cm -s 

(1 to 10 MeV) I 

i 

(c) f-rays ^ 

577 MeV/cm^-s 



1 MeV 



RADIUS (cm) 

FIG, 6-4(1), NEUTRON FLUX DISTRIBUTION 
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Particle I Fl\ax for 1 mr em/hr close | 

I L 


2 

MeV 

1 

MeV/cm^-s 

4 

MeV 

1 841 

2 

MeV/cm -s 

6 

MeV 

1 

2 

MeV/cm -s 

8 

MeV 

j! 

J 1022 

2 

MeV/cm -s 


using these the total aose rates at different layers 

have been calculated & listed in the CODE results. 

6.3.5 WEIGHTS OF SHIELDING MATCRIADS 

In a marine reactor the shielaing should weigh 
minimum possible while still being effective^ once the shield 
dimensions have been arrived at after having analysed many 
combinations of varying thicknesses in the CODE, the approximate 
shielding weight has been calculated using simple geometry ; 

techniques. Por’-the varying heights of different shield -layers, 
varying height- scaling factors have been arrived at, after 
extensive study of a comparative nat\ire done on various marine I 

i 

reactors (information of only tho-se ' available J. i 

The overall shielding weight (including water- | 

snields ) comes to approximately 838,2 tons. The weights of the i 

I 

primary system, Core, Control Rod drive Mechanism (CRDM) and coolant &; 

i 

piping are to be separately considered for the total reactor | 
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weight. An uncertainty error of 10% in the above determined 
value is anticipated' ■ . 

6.4 COMMENTS 

It has been demonstrated from the CODE output that 

the maximum dose rate from the secondary snield is q ^735 mrem/hr, 

3 

and that from the primary shield 4.331 x lOmrem/hr. Overall 
shield outer diameter is 30.62 Ft and weight 838.2 tons/ a 
200 cms space for the primary loop machinery has been assumed. 

Witn an increase in this space, the shield weight will increase 
substantially as the quantity of lead required increases. 

The neutron flux at the core edge has been chosen in 
a conservative manner to ensure proper shield design under all 
operating conditions. Fluxes obtained from chapter-5 have been i 
incremented by 2.304 times (1.269 Engg. hot-channel factor; , 
1.178 control-rod peaking factor- Z Sin(Z) over Sin(z); 1.5408 , 

flux increase due to core burnup) . ^ 

The use of point-kernel technique with build-up 

factors, however, does not fully account the fact that in-scatterec^ 

V-rays is of a lower energy than incident '( -rays, E^.. The scattered 

component consists of a spectrum of energies from O to Ey. Tni^ 

changes with distance traversed through shield. Further, the actual 

scattering process is anisotropic. So the dose rates obtained by 
slightly 

our analysis are/hlgher than the actual rates. As f-ray att€nuaf:ion 

is correctly described by the Boltzmann's transport equations, 
improved analysis results when transport codes are used. 
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In thc' of a mar in- reactor/ the control 

ana instrumentation system is a k^y factor in pro'/ioing an 
efficient safe ana reliable plant. warsnip nas stringent 
manoeuvring reqXiiremonts wbich imposes large load variations in 
relatively Siiorter times. Consequently, tne system snoula be 
stable unaer all conaitions of operation with no unallowable 
overshoots, jvnowleuge of core reactivity balance ana temper atu 
coefficients is reguirea to assess tne above. 

The ixeactor Physics CODb and subroutine AKhOi'I 
in specific (p.l6 of coinj_uter print-outs) have been referred to 
throughout this chapter. 

7*1 Kh. ACTIVITY cpx iPUTi 111 ox 'S 

Two iiUjjortant fission products Xenon & Samarium 
absorb tixermal neutrons substantially. In aauition the long-live( 
fission products also absorb thermal neutrons. Thus, if the read 
were to remain critical it siiould have sufficient reactivity to i 
override tne a.cove deleterious effects. In general, the only : 
significant effects these poisor^have on is a reauction of 

' . . . I 

tficruial utilization, f to a lower value f . Tne poison reactivit 
, is 

- = -ilp ‘ f 


. . . . Eqn. (7.1) 
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The variation of above parameters with temperature 


Siioula be con&iac;rea. 


135 149 © 

7.1.1 ke ana Sm Poisoning : 


Tais is the most iiuportant fission proouct having 
a neutron aosorption cross-scrctiun of 2*7 Hb (2200 m/s). Tiie aecay 
cnain for tais isotope is := 


135 

Te 

/S 


< 30; 


_J3 

6,7hr 




Direct Pist-ion yiela=0.061 Direct Fission 


p 2 , 63 s yr 


( Yj) 


Yiel2±a0^003 

<v • 


(stable) 


Detailea explanations of differential equations involving the 
above caain are given in many boors ic v/ill not be attempted here. 
The following results, nowevor, ensure;- 


(Yj 

K * 


£qn (7.2) 


wrier e, ^ temperature depenaent function 

= Average fuel fit.sion cross-section 


= Average core thermal flux 


= Steady-state Xenon Concentration 


From Ref (23) 
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^ 135 

Eqn. (7,2) substituted in Eqn. (7.1) (forZT^p) gives the Xe 
activity at steady state. After reactor shutdown, Xe activity 
rises above this value. The maximum occurs at a time, t 


max 


max 




In 


( 


A 


X 


1+(Y^ + 

Eqn. (7,3) 


, =X 

I I 


where, /a 


ax 


The Xe concentration is tnen given by substituting ^ 

following Post- shutdown equation 


t 

x(t) = Xq e 




Xj - 


( o 




<^1 + % 


where, ^ g 


ax • % 


t 

e ^ ) .. Ecpi, (7,4) 


, Steady- state Xe 
concn. 




= 




, Steady-state I 


135 


concn 


obtaining from Eqn, (7.4) and substituting it in Egfi. (7,1) 

3X 

the maximim Xe activity is evaluated. 

149 


(iD-saa- 


"Ihough a lesser poison, this has a thermal absorption 


cross-section of 40.8 Kb. It *s decay chain is 

^ 


2 ,hr 




54 hr 


149 

Sm ( stable ) 


direct fission yield = 0,0113 

«p ) 
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The Steody-state Sri absorption cross-section, * 

is ylven by 




1^9 

Unlbce Xe txiis is inaepenaent of neutron flux. Tne max. Sm 
concontratiun after snutaown is given by 


S 

max 


S + 
a 


P 

o 


bqn (7*6) 


occur iny as t- — 5^inf, ’where 

^o~ ^f ^ aS i subscrijjt ’p' stcands 

- 3! 

P^= ^p^f ^ Promctnium. 

Obtaining from i;.qn,(7.6) anu substituting it in ii.^.p(7.1)/ 

the maximura Sm activity is evaluated. 

7.1.2 LGhG-LiysP FISSION PuOOUCTS 

Tnese also absorb tuermal neutrons. i’'iany loroaucts 

naving varying cross secti>^ns for neutron absorption ana decay- 

25 

'times ensue by U burnup- (Refer Table 9.2,Kef (24) )j; for 
computing tne long-term effects of trie accumulated proaucts, 
a-n average cross-section of 'iS b per fission is chosen 

(p. lObykef (24) ) . Tiiis has been obtainea neglecting the rapidly 

’ V- 

satTnrating proaucts.- Since the. cuncentrati.. n of ‘tiiase products' 
i§ essentially tiie bxirnt atoms of it’s macrpscopic crotos- 
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section/ is given by 

^afp = ^frac * """ * ^<J"- <7-’) 

25 

wfiere/ Fraction of U burnt at End of core life (see 7.1.4) 

n^^ = nui'iiber aen&ity of U'^^z :^/cm^-s 

Substituting Eqn (7,7 in Eqn. (7.1)/ the activity of the above 
proQUcts is ootainea. 

7.1.3 Tu:.rip^i:s.^.Taioi. CCEr j'lCIEi.TS Qi?' i<_.fiCTI-/ITY 

Tne temper atiire cocsf ticic-nts of roautivity ( '^;p) are the basic 
parameters rev-^uirea ror xveactor Stability i-nalysis, 7iS 
changes v/itn teiuperatures of operation/ a knovvledge of ' 

s 

at various reactor operating Ivels both for tiie fuel & moaerator 
is essential. Unlive land-basea reactors/ tne load variativ^ns 
here are trora full-load to below 20% for automatic control. 
Cvjnseg.aently/ ^^'s at various loaas woula be difxerent and require 
separate calculations. For convenience only full load, 3/4, h 
ana ^ loads, as well as full-load at EOL (Ena of Life) have 
been considered. 

The basic equation for°^;j, is 

+“^1 ‘h > 

^y(i) = 1__ ^ di_ 

i al? 


Where , 


is^ru (7*9) 
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i = Paramc;ter iuvulved 

~ Tnernial Non-"lea,;ac;e probab±i.ities. 


Eq.n(7.cs) can be expanded as be low: - 

=._ T . (b ) { 3 ^) ) 

1+3 

+C\^(3^) ) 


Eqn. (7.10) 
Es^n (7.11) 
Eqn (7.12) 


Eacn. of tile terms in Eqn. (7.10) to (7.12) can furtner be 
exi:/anaed but not aone here. (Rarer Cnap. 13 ,Kef (23)). 

Trie basic assumption for coiiputations here is 
teiiiperatures calculatea are uniform for tne materials’ involved/ 
tiirouqhout the core. In tne Reactor Pnysics CCDE, every parameter 
is corrected for the operating temperature (except tho core 
airnension variation wnicn is negligibly small), iience is 
aeterininea by giving a small perturbation to tne operatiny 
teriiperature ana aeteruiining the consequent variation in 
and nence^Sp using Eqn (7,9). Tnis was aone for temperature 
perturbations of botn fuel and moaerator operating -temperatures 
inuividually , giving*^ and tJrie respective tempercture 

coerticients. The average fuel ana moaerator temperatures of 
different loads were obtained from the Tnermal Hydraulics CODE 
(dnap 4) by varying the linear power rating of fuel rod 
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api^ropr lately. Results obtained are summed up in Taole 7.1. 


Taolfc! 7.1 Reactivrity _ Coefficients of Fuel & Moaerator. 


1 LoaD I 

^frac 

1 

jv ^ FUii» 

, TEi'^PK 

1 

Full ! 

1 

; 

0.b>197E-3 ! 

1 

0 • b>929^ 

508.07 

3/4 1 

i 


0.yy2v5 

442.27 

1/2 1 
1 

l.lt565Ji»«5 1 

j 

0.ys2&6 


1 

1/4 i 

i 

, 1 • 32l4£*»3 j 

1 

0. s9i08 

1 1 

329.43 

f 

1 

jci/OL, ! 
Flji-iXj ' 

0.0 

1 

! 0.sy2y4 ' 

i 1 

1 J 

508.07 


.^V. MOD 

j TEaPR 1 ^ 

, . 

; 267. 72 

{ 

1 

2.7376E-4 

i 6.23273-4 

1 

i 2b5 ^ 84 
i 

2.8963b~4 

j 6*04t64£— 4 

1 

j2d3, 93 

1 

3.C667E-4 

i S.6/J3E— 4 
» 

i281.y8 

i 

3 . 24683—4 

1 

1 5.320e-4 

■267.72 

2. 3940 E-4 

i 

; 7.00E-4 1 



Z.*. r Burt j up^___ Sc R e activity 3 alc^c e, 

Knovvleuye of tixe fuel fraction biurnt is necessary 

to obtain tiie core burnup and reactivity balance. Tiiis nas been 

acaievea usii:iQ a term ( '3UKJMT ' in CODE), ie that fraction 

25 

of burnt U ^ below Wiiich the core is ui.^ble to overriae tne 

Xe-aii ana fission proauct activities and thus becomes subcritical 

with all coiitrol roas reinoved (P^ =: 0.0 ). Tnis state of core 

frac 

is termed EOL (End of Life). 

Tne core inouilt reactivity, reactivities required 
to come up from cold to iiot-zero- jpower^ ana zero to full- power 
are all aeterinineu by ciioosirig tne appropiate temperatures of 





fuel cc moaerator in tne CCDn: and evaluating .c . «t the t.cL , 
core burnu^_' is given by the relation;- 


Burnup . (wt. of U per ton ot UO^ * 3^^^^ ) 

(i'i>Aa/ton) , , , , ,,i5 , ^ ^ 

Vvt, of U ournt for producing 1 Maq 


. . bgn ( 7 » 13 ) 


ana, finally, txie core endurance (in days of full power ) is:- 


rlnaurance 


Burnup^^*^ JTptpl wt . _.pf„up. 2 .. 
operating power level 


. . . . £;qn{7.1'y:) 


(x<e±or line nos.69CMD to 7400 o£ p^l7 of computer printouts for 
i!,gns(7.13 ) and (7*14) ). 


Having known tne aleove results, tne reactivity 

balance of tne core is rnaae as in Table 7,2. Tne aiaount of HoISOiSi 

required in tiie control rods snoula be suen that tne core remains 

subdr itical even in the Cola condition. Honce, a P =2.35i!;-3 

frac 

was found adequate to drive the core subcritical by over 4% 

Ak (:c^^^=0. 9532) - Tnis means as tnere are 6 control rods per 
assembly, by a sim^-le geometric coifiputatien, the rod aiameter^ 
is 0.51 cm (J.494 cm 3^C pellet aiameter), assuming a stainless 
steel claaaing of tnicjcness 0.53 mm • 
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7. 2 CORJ:.. i 

Cola to 'O’ jjO'vvtr 9.6360 

*0’ to full 2 -:.ower 5.4634 

I jie, Sm poisoning 2.8563 

i 

T^O/Sm ovorriat(at EGL) 0.7356 

Fission Froaucs Jr-oisoning 1.8577 

Burnup coiisiaerativ. ns 10.2373 

Total iiibuilt reacti\;-ity .... 30.7863 

i 

{ Total Control rod. v^ortn 34.93 

Burn up (/iVorago) - 146-57..i8 MWd/ton 

£.naurance (days of full p®wer) - daya^i 

! 

Till nov^, tiie basic reactivity reguireutentS/ 
tfaiiij^erature coeiriciants 6c core burnup nave been computea. 
Next tne control systeuii will be evolvea ana later certain 
stability analysis done. 
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7.2 _SYST£/i 

7.2.1 S a lient Features^.’^-:. 

Tile,- raanoeuvrin^ of warsxiij^s iiUj_oots larsae 
load vexiations on tno reactor, ii conveiitional control is 

inaac--‘_uatc; for ti.is purpose as it is aesignoo for only a small 
variation from full loau(Ui::.to lO/o)- tno basic signal it uscS 
is tiiO slow resj^onoxi.g avera>^fc coolant temperature cxiange. In 
marine -reactors/ nowe'/e-r/ an aooiti. nal j^roiapt sivynal will be 
required (rief (32) ) v;nicn 'feels' tbe i-ower trarisient as soon 
as it startS/& aisaj-.'pears graaually at tne ena of loaa variation 
after w.-icn only tbe slov/er reSx-onding signal would oe active. 
Care saoulo be ta-.ei that the reactivity insertion rates ao not 
maiN.e tiie reactor proaipt critical. For tais, the rate is restrictei 
to 0.03 Vo A 3'^ (i.o. axj^xoit. 0.04o^) based on the fact that reactol 

snould respond by a 2.3Vo/s x/Ower variation to go from S<0Vo to 20% 
in 30 Sees (ko±(32)). 

TiiuS/ in tile first part of tne ox-eratioual ; 
tra.asient/ the control system siiould act promptly aiia efficientli 
performing gross regulation and later fine regulation to prevent 

f 

it ax-x-xoacnin.;^ unsafe conaitions ana oscillations ai'o^na ■ 

egualibrium. 1 

7.2.2 Tn e^ SX,sJa=-£n : - | 

1 

i 

1 

I 
I 

I 

! 


A simplified aiagram of the system is snown in 



POWER CONTROL 
SYSTEM 



RIACTOII CONiabL SYSTEM. 










Fig. 7,1. It consists of 2 subsystems/ viz. 'Po^jer control systerri* 


(prompt acting) and 'tomporoture control s;ysttra' (slow acting), 
(i) Poy.ar _Co/.brol Systoins- 


nerc tne main signal is inbalanct. between 
seconaary steaifi flow rate, rh / and core g>=n6ratea power (proportional 

o 

to neutron flux, 0^,) . 2 aaaitior.al teruiS/ d>-.rivcitives of tne 

aforesaia, are aaoed for enhancing prompt cntrol, v/itn scaling 

factors 0. cc fi . Trie rt^sulciny -over error signal activates 
i a 

eitiior control rod cluster u-idC) x, B or C aepenain^ on tile signal 
iiiayijitUQe. -■ is usaQ to sanse transient coiioitien; if signal 

Qu 

exceeds a CorLain val ue_, D3 (Dead Band) is cuanped from 

(wiue/ normally) to SP^ (narrow). Rod speea is coatroilea by an 

erior signal siiailar to tiie above but witii a different gain factor 


b * The speea Control unit assigns a minimum & maxii’auifi insertion 
rate ( = 2.3 >o/s j-KJWer variation) as sixown in the fiyure. 


(ii) Teiapor aturp_ CpntroJ.^_ Bp steia i 

Here maiii signal derives froai aixforence between the _ 

primary coolant average teaii-r. (T, ;) and a constant reference tea'ipr, i 

av 1 

1 

Asoent for rod~iaotioii is given wuonever error signal e.>^ceeas tiie ’ 

I 

tempr.DB. Trie cnaroct eristics of tiie various units are as siiown | 

in Fiy. 7 . 1 . I 

[ 

I 

Signals o>/tained froai tne' above 2 su.- .-systems 
are suamea up to obtain tno nett error signal. Trie sequential 
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co.itrol unit afecidts tne kcC tnat is to be moved, \vith tiie speed 
assigned by tlie rod-si>eeQ control units. Tne control-rod actuators 
activate tne- particular ROC, tiius finally varying 0 ^; tais feedbacrc 
is reflectLu in tne pouer control system resulting in a reduction 
of tne original error signal; tais effect continues ensuing 
automatic control. 

7,2.3 Control Sen erne 

i-it:re, a brief coijcept of tiie roa-coatrol 
strategy ov^r tn>^ core operating lire time is eluciaeted. 
ii) Control group is constituted of 3 RCC’s A, B ano. C as in 
Fig, 3.3; eacn cluster moves as a bundle, inaependent of eacn 
otaer . 

(ii) A constitutes tiie 'Sciam’ rod, capeible of sriuttin^ down the 
reactor ana is jvept fully retracted unoer nonnal operation. 

(iii) B constitutes 'coarse' and C, 'fine' runtrcl clusters, 

(iv) Initially, B is Jeept fully inserted ana C moved till it 
reocnes tne top limit. 

(v) Tnen C is fully insertoa and B half extracted; C is moved 
up till it j-eaches tne top limit. 

(vi) Finally/ C 'is fully ins>_rted & B fully extracted; C is moved 
till it reacnes th^ top limit ( Txiis tvoula taen signify tae end 

of core- life). 

Tne above scxieme is only a concept & saould 
be supported by extensive calculatioxis wiilca is beyoad tne scope 
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of tnc prc,sfeiTfe Vv-oxk. 


7. 3 STi-^lJulTY Ai'./.LYoIS 

7.3.1 Cpmputatiyii^ 9± _somc; parami:st5;rsj__" 


(a) Tne averavjo neutron life-time/ 1* is 


2 2 

1+ L ^ B 

1 *• .t- A — 


V 


Eqn (7.15) 


^vHero , V 




a 


= neutron spov,a (average) 

= Fuel aj^sorption cross-soc (avera-^fe) 


1* ciianges witn load as 2* cnanges vjith operating fuel 

3 . 

tamper atures. Taole 2.3 gives 1* for aifferent loads. 


(b) For neutronics calculations, one-~precurs- or group-moael 


nas been cuosen. Txie transter function, is given by j- 




r, 


n ( s ) /n, ^ 


s+ 7\ 


Egn. (7.15) 


<r (s) 


lax (s+X+p/l') 


Wiicre, X= 0.013 4" 
p = 0.0065 
Rest of tiie 


one-group 
J values . 

symbols have t.^eir usual meaning. 


7.3.2 Two- teiCipe r ature Feedback : - 

Txie control loop ror toaiperature feeobacfv 
effects of fuel and moaerator is 
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Kn ) 


( X ^. S ^ Ca ) }s — 


a„.Gr (s) 
m m' I 


6ii(s)/ii. 



Wt; JTfciCiUiirc.' to cVcilustte tilO iflOCitslT StOlT , flifel tlTcJ.nsf'-:^ functioiiS/ 

G ana 
m £ 

Taking an eloinontary nuat balance for a fuel rod in tne reactor 
anu a^^i-lying a -small perturbation tc tne fuel temper atxiro/ T^. 
We arrive at the tollowin^, relation (from p.92/ ref (34)). 


bgn (7.17) 


1+S'i^ 


For tne iiK^aerator a similar e^^Ueti. n as aoove is obtainea:- 


bqn. (7.1b) 


1+ s^m 


where, At^ = (T^-T^) at steady -state 

/}\T = (T ~T_) at steady state 

m m C 


= fuel time constant = m^.C^^/h^ 
'h£= Pov/e-r produced per fuel rod 

■ A * 

At^ 

t = iTKiaerator time constant. 


Refer AppenaxA a7.1 for evaluation of ^ ana t 


•..*v C 
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Table 7,3 Two Temperatxare Feedback Parameters 


Full 

Load 


5 ' 


3.693 i 3.409 


3.162 


2.963 


3.693 


> o 

At^ ( G) 


220.19 156.44 99,05 


47.46 


220.19 




0.483 


0.482 


0.481 


0.480 


0.483 


At ( sec ) 
m 


7.72 


5.84 


3.93 


1.98 


7.72 


H (0) 


0.0651 0.0488 0*0326 0.0165 0.0581 


0.0263 0.0206 0.0142 0.0074 0.0255. 


1.3882 1.4421 1.5088 1.5664 1.2788 


1 (ps ) 11,31 


10.55 


9.82 


9.11 


1-7.43 


Points of Intersection on Imaginary axis of G(s).H{s):- 

w (H^) + 27.74 + 28-02 +28.31 + 28.56 + 22.33 

GH (xlO®*) - 6.252 t 6.169 | - 6. 126 ( -r 6.073 | t 11.956 

+ + 
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FIO.7.2 NVQUIST PLOT OF 2-TEMPERATURE FEEDBACK SYSTEM 
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7.3.3 Boiler Transfer Function 



Tills figure snows 
the lumped moat;! used. 


Assumptions niade in 


( i ) Pr imar y Temper at ur e 


% ( T^,-+ 


^ ' "bi bo' 

(ii) Secondary Temperature is that of saturated steam, 
wnich is a function of pressure. 

(iii) Steam-water mass in secondary side is constant. 

(iv) Throttle is kept constant; power output of boiler is 
proportional to secondary pressure, p . 

using basic heat balance for Boiler Power (BP) the following 
result ensures (from p.99, ref (34)) for Boiler Transfer Function 




Eqn(7.21) 


i€iere K = 
' X = 


1+ ^+S 
X— 1 


X+1 

2 ' ^ 


V P; 
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, <ra C ) secondary /BP 

T;+ ^ V- 


+ »./P 


s' -^S 


^ == 


(mC ) secondary / BP 


In the above equation p is the slope of saturation curve , 

3 

(d p /dT ) ' at the appropriate pressure of operation; it is got 
s 

by purturbing the following equation (p in p.s.i. ) 


T (p) = 235.7 + 0.061294*p - 0.5661E-5*p . 

s 

The secondary side pressure has been chosen as 
500 psi at full Idad & proportionally lesser for part loads thereof 
(by virtue of assumption (IV) above). As before boiler transfer 
functions have to be evolved at various loads and results are 
in Table 7,4* and 'C are still unknown (as secondary side 

design has not been due). The Bode plots, however, for Eqn. (7.21) 

may be plotted, but break frequencies cannot be determined. 

* . _ . , , . , , . T, -r" ✓ j.' which is a design 

As w — > xnf, the plots tend to K~^ /t;+ 

parameter^ and at low frequencies, the plots tend to K. 

Fig. 7.3 gives the plots obtained. 

Table 7.4 Boiler Trartfer function parameters . 



Full load 

75% 

50% 

25% 

T 

S 

240.688 

255.732 

204.977 

173.524 

pu/^ (°C) 

S S 

71.256 

53.255 

48.829 

42.109 

K 

0.7675 

0.9021 

0,9404 

0.9741 

K. t ” 

% + i 

0.628 

0.823 

0.905 

0.964 
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Fiy.7.3 Boiltjr Transto-r Function” Bpd^_ Flot__ 

7.3 _Rjb;^SULT_S _ , COid'ipi'JTS 

”Tne eofc R>_oCtivity Balanct; is aisployed in 

Table 7.2. From figs. 7.2 & 7.3, tno reactor is stable unaer 

all conditions of operation. Many parameters are xiowcver still, 

unicnown a re-.iuire extensive calculations to be deterirtined- In 

tne reactor control system and all dead-bands are 

yet to be evaluated & lienee an overall study of trie system would 

not be xjossible. Transient Analysis & reactor dynamics iiave 

not been done in the present study but do require careful 

consideration for erfiedk of unallowable overshoots. Throttle 

variations as described in Cnap.4, ref (34), were studied 

to an extent? but further progress was not possible due to many 

unicnown parameters beyond the scope of the present design. 


'k'k'kic'k'ki^'k'k'k'k'k'k'k'kiK'kic-k'k'^ 



CHAPTER -8 


RESUME 


REACTOR ON SHIP 


8.1 SHIP'S MOTION EFFECTS 

8.2 COLLISION BARRIER & GROUNDING PROTECTION 

8.3 MISCELLANEOUS FEATURES 
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Inte'i^rt-tion of reactor siilp 

includes Certain specific doaipn con siav_-riit ions xiitiiorto not 
discussed liice snip's motion erfuctS/ induced vibratioxis, collision 
ci grounuiin^ protection for reactor^ etc. 

0.1 Si'ilj: 'S i^pTlv-’is Jii; Oi'i Sla^'i iroRFOidlidsiCn^j — 

location of tde rcoctor on bocird snoula be sued 
as to mii'xiinizfe snip's roll & piten on the steam generator- ' 

location at about tno snip's CG(Contre of Gravity) would be ideal. J 

In the U-txxbo stoani generator ^ faea watar gains 
boat by a 2-pnasu natural circulation system- firstly/ feea water 
absorbs tnermal energy from primary coolant becoming partly steam; I 
nex-t/ tnis steam-water mixtxxre from the t^e l*undle passes txirot^ [ 

a steam swirl vane locateo in the riser whicii separates water i 

i 

I 

from ste*am cert if ug ally and subsequently through adaitional 
dr yers to decrease the moisture content. Txiis 2 phase system/ 

I 

under the effecr of snips 3-aiiaensi./nal motion could present 
flow oscillations causing loss of heat transfer efficiency/ f 

large level fluctuations in the down comer/ pulsation in Steam 
Supply/ etc./ thus forcing a barely stable reactor control 
unit into a region of instcibility. Special siiiulation coirputer 
CODiiS have to be developed to study the etfects of snip iriotion 
on steam generator circulation; tnis is not attempted here but 
a brief sketcii of the associated problem given. 

1 
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The recirculation loop is divided into a numbt^r 
of small equal portions of, say, Az, The normal heat balance is 
done for each element, Additionally/ the rolling & pitching effects 
sxjbject Az to tangential & ctntrifugal forces/ modifying the 
momentim equation. The dipjAng effect results in a modification of 
'g' (accln. due to gravity) by an oscillatory motion. Consequently/ 
the thermodynamic stability of the recirculation loop in .presence 
of the above effects is done. Finally comes the secondary water 
level control system. Power level/ feed water flow & secendary 
pressure dynamically alter this level; tfxis necessitates the 
development of an auto level control loop to ensure system 
stability. 

Extensive studies in the ship motion effects 
have been done by iSUSocolo, G. Selvaggi et al (Ref (37)) using 3 
main computer CODES- Sabre-3M, F-111 MO and DCR-2. The study 
evaluates dynamic behaviour for following perturbations:— 


Roll 

30* (max) 

; period 

= 12 s. 

Pitch - 

10* (max) 

; period 

• 

CO 

11 

Dip 

0. 3g 

; period 

= 7s, 


Results obtained reflect the overall performance of a marine steam 
generator as summed up below:— 

# 

(1) System presents an oscillatory performance for an 
impressed step perturbation, . 

(ii) Insertion of an orificed plate suitably positioned in 
the downccMner region increases stability margin & damps the 
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oscillativ^ns Inducea by snip's motion offsets, 

(iii) Strain <_j,uality incr^asos '..itn- docroosini^ wator IcVol in tno 
doVi/ncomt^r . 

(iv) Systoin stability aocroas^_s v/itii incroasinp stv^ara gonorator 
inclination or lowor system op^-rating prossuros. 

(v) SystoM stability incr^-c'-Sv^s witn lowor pr._ssure crops in noatod 
zona or lowor £o_q w<.to-r tomp^^rcturas. 

2 CG i -iXjI S lOn ai'JL) GkOo; jbXiJG L-'ROi'rjCT lo X‘'l 

To miniiiiizt collision otfocts for reactor, 
protection is roguired to avoid ro._ctor aamago; conseguenting in 
uncuntrollcQ roloas<_ of fission prooucts (approx. 6 ciirios/V\/ 
g^-noratod ) . T..is is aone by a grid structure of about 80 cm x 
yO cm arranged betwov-n longitudnal bulicdoads & outerhull/ 
supi.ortfed by decks benind the longitudnal bulkao-ads Cref(3o)), 

Tnis structure transfers the collision load into tile transvoxse 
bulkbeeas of tno reactor corapartment , tnus destroying tne stri-ving 
blow. As per a test progrein conaucted (p. 422, ref ibia), no 
mecnonical damage will occur for components benind tne collision 
structure; in addition, the protection is adeguate against 
penetration of a jet aircraft at nominal speed. 

For the grounaing protection, a grid type 
support structure 140 cm i:iiQn is located above a 200 cm niyn 
double-bottom (DS). If db is damagea, the grid is capable of 
witnstanding the weight of reactor containment. Conseguently , 
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during severe loads# the DB being more easily deformable# acts 
as an elastic spring. 

8.2.3 MISCELLANEOUS FEATURES ; - 

Vibrations are excited by sea-motions# 
oscillatory propeller forces & axixiliaries; frequency range varies 
from < 1 Hz to 80 Hz. High energy vibrations from propellers are 
specially significant & require detailed analysis. As compared to 
conventional ships of equal power# the machinery weight as well as j 
the normal cruising load on machinery is nigher by 1.7 & 3.5 times i 

respectively (approximate figures). Trie Power consumed by the ! 

! 

primary coolant pump is substantial. | 

The ship's speed depends on the SHP( Shaft 
Horse Power) developed by the engine. The speed versus SHP 
characteristics of a ship is dependent on the displacement# profile 
at the water line & certain other factors; it however follows an 
approximate cubic law with linear speed increments requiring cubic 
aiP increments. Consequently# at 50% of power developed# the 
ship acnieves nearly 75% of its maximxim speed* A rough estimate of- 
22 knots may be expected as the maximum speed of the daip presently 
under consideration# with a 30 #000 aiP maximxim power output ( assuroinc 
70% of steam to be used for propulsion only^ and a plant efficiency 
of 22%>. 

' An important aspect hitherto ignored is 

requirement of a secondary propulsion system. Coupled to the main 
gear box# a compact gas turbine of approximately 15000 SHP may be 
included so that# in the event of a main engine breakdown# the ship 
is still capable of making way at 65% of its maximum speed which# 
now/ would correspond to that of 45000 SHP (*^25 knots). 



CHAPTER -9 


CONCLUSION 


RESUME 


9,1 SUMMARY 


9.2 SUGGESTIONS FOR FURTHER STUDY 
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By now, the salient aspects of marine 
reactor design have been covered. In conclusion, the destroyer 
class ship is about the smallest to wnicin nuclear propulsion can be 
reasonably applied - a preliminary study of containment weignt revea 
this (chap 6), as a propulsion machinery weight of ri; 2500 tons is 
bound to accrue ( & a displacement of the ship of 8000 tons). 

9.1 SUMMARY 

Some of the design parameters evolveO are ' 

t 

briefly summed up in Table 9.1* Sinct^ tne vari^w^us design processes ! 

f 

are mherently interactive, change of a parameter may nave a | 

retrospective effect in other processes. In the shielding design [ 

I 

a macninery space of 300 cms (instead of the 200cms assumed) would j 
have resulted in a dose rate ofo,4982 inrem/hr outside the 
containment with a weight increase of 262.0 tons and diameteir 
of 6»6 ft. This space can be fixed only after all components of 
primary/secondary loop are designed. Core burnup obtained has been 
low for a 5% enrichment. Tnis-immaiation is due to the 

compact, close— latticed core configuration cnosen that is so 
vital in reducing overall containment weight in marine reactors. 

Table 9.1 Design Characteristics 

I GENKRAD REaCTQR DATA 

A. Reactor Thermal Power .... 120 MW 

Mo. of Primary coolant loops ... 2 
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3. Primary Coolant 

- System prfcss\ire 2000 psi 

- Total flow rate 5553 tons/hr 

- Inlet tempr, 280 °C 

- AVer ago rise in core ...... 15.42 °C 

C. Fuel (VO 2 ) 

- No. of cere regions 1 

- No. of fuel asslys. 29 

- No. of fuel rods per assly. — 282 

- Enricnment 5 v; % 

— Core Power density .... 63.51 MW/m 

- Eqbm. core average burnup ... 14^697.8 MWd/t 

- Wt, of UO^ In core .... 6.919 tons 

D. Control rods 

m 

— Type/Niomber RCC/29 

- Absorbing Material B^C 

- Rods outor dia/pellet dia ... 0:< 51/0 ,494 cm 


E. Secondary System 
- Steam pressvure (at full load) 

T Steam quality .... 

F, Plant Physical Parameters 

— Overall Wt. of shielding 
— Containment outer, dia 

ii» nuci£AlR design data 

A, Structural characteristics 

- Core dia/hei^t .... 134.29 cm/133. 4cm 

- Fuel rod & control rod cladding .... ss347 

- Fuel cell type/pitch ... square/1. 30cm 

- Fuel rod dia/pellet dia ... 0.99/0.858 cm 


500 psi 
0.9975 


838.2 tons 
30.67 ft. 
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-Voliime fraction H 20 /UO 2 in oort; ..... 1.592 

B. Piiysical & performance cnar act eristics 

- Fuel type/density sintered UO 2 /IO. 97gm/cc 

- Fuel Enrichment - 5 w% 

- Hot channel factors j 

F^-t / F^ /total F. 3.44/1.04/3.58 

Q Q w 

fJ / / total F „ 1.994/1.22/2.105 

R li -rt 

- Flux functions analysed z-Sin (z)/Sin(z)/(it-z)* 

DOL/Normal/EOL Sin(ix-z) | 

C. Control characteristics ' 

- Type mecnanical, RCC j 

- Absorbing material .... B^C(den 0.7 gm/cc) 

— Rods per assembly /Ne .of . asslys. ... 6/29 

D. Kinetic characteristics 

— Moderator tempr. coeff. (full load) . . .- 6 . 233E-4 

Ak/k 

-Fuel tempr. coeff. (full load) . . . . -2. 738E-4 ^k/k 

III. Core Reactivity Balance 
Refer Table 7. 2 . 

9. 2 Suggestions for further study 

Some important aspects of marine reactor 
design not done in the present work are;- 

(a) Materials design - The performance of only a few materials 

has been studied; optimization study of different reactor materials 
is suggested to evolve better design parameters. 

(b) Ship-motion effects on the steam generator require a study 
in far greater depth than mentioned in chapter 8. This proolem 
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is peculiar of marine reactors and has a direct effect on overall 
system stability, 

(c) In-core instrumentation Cesign requires careful consiaeration 
as it conveys core fluxes and temperatures at various points in 
the core; these signals structure the basic inputs for the reactor 
auto-control loop. Control rod worth evaluation is also suggested. 

(d) In the Snielainy Design tlie radiati^^n from reactor effluents 
requires evaluation after steam-generator and piping design is 
done; its magnitude is, however, siaall & tne shielding thicknesses , 
may require marginal modificatien to counteract this additional dosf 

(e) The pressure-vessel & piping design for mechanical, thermal | 

and other stresses necessitates careful checking out; in the | 

present design a vessel thickness of 12 cm has been assumed I 

I 

based on available literature and not on calculations. j 

(f) Design of press'urizer , steam generator, components of | 

I 

secondary loop & subsequent heat balance is also proposed. 

(g) Reactor start-up & shutdown procedures are suggested for 
the present design. 

Study of an extensive topic as a reactor 
design In a limited time frame effects in simplifications & 
exclusions. Ihe core flux functions for thermal hydraulic study 
do not precisely represent actual flux variations. Other flux 
fvinctions as considered adequate may be assumed to evolve more 
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realistic performances, llie core nuclear design (Reactor Physics 
& Shielding ) has been done using a heterogeneous 2-group theory 
with cv-nstant core temperatures for fuel & moderator- 
multigroup theory using local temperature variations and 
development of transport theory computer CODES is suggested; 
this evolves a more precise design/ closer to a practical 
reactor unit. The reactor automatic control loop may be studied 
in detail/ both for system stability & transient characteristics/ 
once the loop constants (of Fig. 7.1) are evaluated. Study 
of reactor transient for an ahead-astern simulation consisting 
of a load variation from 90% to 20% in 30 secs and immediate 
return to 90% in another 30 secs would be of significance in 
ascertaining the system performance. 

An earnest attempt has been made to describe 
the complex# varied & multi-dimensional problem that a warship 
reactor design is. Both strategically and tactically# nuclear 
propulsion is the ultimate form of propulsion for warships 
today; it is obvious to all who pause to ponder this matter. 

It will remain so lintil a superior form of ocean propulsion 
for suppcjTtir^ modern weapon systems is architectured, 

i^ikicicicicicic'kic'^icic'kicic'kic'k’kicii'kiki^ 
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— (zSinz) =0 at maximum ; 
da 

■=> Tan(z) = -z 

The above eqn. is satisfied when z=2, 02877 
So/ maximum = 2.02877*Sin- (2.02877) = 1.8197 


We shall choose the axis as above in diagram. 
H 


^av 


H 


r z. SI 

U 1.6 


Sin.z.dz 


where, e= 


H -H 
e 


8197 


As z is distance in cms we will have to convert it to radians* 


= % 


H 


^ H 


^ * Sin( — ^ )*dz * — i 


H 


e e 
7E „ 


H 


^dx 


% 


H 


.dz 


1.8197 


Let 


X 



15b 


1.8197 * H 


X. Sin x.dx 


= ^ 

1.8197 


(,- ^ ) i- 


; HU ^ (.HU ) + 

i - h7 — 


-. Cos (- 


•) + Sind 


Sin (- 


Putting A = 


B ~ ^ 
' ® ’ H 


(H -H) U 
e 


£7 


1.8197* ^ * H 

# - 

H^(-ACos(A) + Bqos<B) + Sin(A)-Sin(B) ) 
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3. TRUiSCATED BESSEL ‘S FUNCT.ON ; 



Function 


F(r)=AJQ(lr) + CI^CnoT) 
R = core radius 


A,C/ 1 & m are all constants determined 

from tne Reactor Physics CODE and hence known. By definition 
H R 

1 




J/ 
0 0 


F(r) * (Axial functipn) *2 rdrdz 


Ihe axial function integration is done in (1) arxl 

( 2 ) above . 

Ttie radial integration is:- 
R 

c^= J ( A-J^Clr) + C- J^Cmr) ) * 2.r.dr 

^ 0 

Since ^ Jlr .S q (I r) .dr = lr,J^(lr) 


and 


s. 

mr*>3^(mr)dr = mr * I^(mr)^ 

>4 


V ^ R.i^(nR) , 


15S 


At r= 0 ^ the peak flux occurs; its value is (A+C) 


'^av 


A + C 


2A- + 2C‘ I^(mR) 


IR 


mR 


1 + 


A 


2 . 


IR 


2( , I^(mR) 


mR 


As 


1 . 


N 


2 J^(IR) 


IR 


^ 2.( f- ).jV^ 


mR 



159 


A 3.2 Di;;TJ:l^IlMlNG 

Let K^^(T) - a+ bT ^ be a linear function of temperature. 

AIM To determine what value of 

(T) is to be usea in the 
fu 

heat conduction equation/ 

Eqn. (4.3). 


Solution 

We shall develop an equation relating and 

using K^^(T) = a+bT then compare this with Eqn»(4,3). 

In the steady-state condition 


q' .Tc-r^ Al = - K^^(T).2'ii;-r.A 1. -g- , in any fuel ring 
of thickness dr at radius r. I ht elating from r=0 to 


^ q ^ ) 

2 Q 

I 1 1 

. . q 

(or) - 


^0 


^ - 2 (aT 

0 


HP 4-'T* 

4{ a+b: ^ O £ j) 


Ihis means K^^(T) in Eqn (4.3) is evaluated at the arithmetic 

T + T 

rttean fuel temperature for determining i.e. at T = _0 f_ 


2 
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A 5.1 calculation OF z- 


Element 

NUM DEW* 

( 

2a* 

^s- - 
(13) 


(l-!io) 

Ztc 

F 

e 

6.044 E22 

0.093180 

11.4 

0.68902 1 

D.9b80 

0.68075 

C 

r 

1.6719E22 

0.030758 

4.2 

0.07022 

D.9871 

0.06931 

N. 

i 

8,2291E21 

0,025058 

17.0 

0.139895 

0.988S 

0.13829 

s. 

1.7536 E21 

0.000178 

2.2 

0.003858 

0.976 

i 

0.00376 


4.1611 E20 

0.000290 

6.5 

0.002705 

; 0.9928 

^ 0,00269 

M 

n 

1,7541 E20 

0.001401 

1.9 

0.000333 

0.9878 

0.00033 

c 

! 

3,2091 E20 

I 0.0000009 

4.7 

0.001500 

0.944U 

0.00142 

i 

Total 

0.150867 

1 

0.90754 


3.896555 


* Taken from TaLle A6,l 




0.318241 cm 







161 


A 5.2 CALOJI^ATION OF D, 


The fast diffusion constant for core, is 


, given by the expression :■ 


Eqn. (1) 


where, ^la ~ av-fast absorption; cross section, 

= averaged scattering cross section, 
since, the scattering cross section of Hydrogen varies 
strongly with neutron energy, we use an artifice to calculate 
its scattering cross-section. Let us make an initial guess 

of ■ iO barns - 

Ihe parameters are 



H. 

0 (total) 

S.S. 

Total 

j 

n ( #/cm^ ) 

0.3641E23 

0.3622E23 

0.8407E22 | 


% <*=> 

_10_ 

CM 

• 

11.4 


'S 

1.0 

0.»)L2 

0.353 



0i333 

0,258 

0.988 


7 . 

s 

0.5094 

j 

0.15212 

0.09584 

0.7577 


Fbr the above parameters. 






t 


0.6862 
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(1 






r.'mc 




X = 

Ws 


ln( §». ) 


S (TOT) 

H 

fck 


= O1.552OS 
% ' 


1 

* tr 


0.7475 


K ) 

We sefe eJDove & one fnorn tabulation do not agtee. 

So we guess hiciaer values of and continue calculations 

We obtain convergence when ^ (H ) = 13,7 

3 ^ 

"^reac ~ 0,69693 




= 0.54459 


= 0.7488 cm“^ Eqn.( 2 ) 


arsac 

Next we have to determine 




N-I 




ln( j 0 _ ) 
^th 


we 


know from calculation for *p ' , that 




= 21,254 barns for UO, 


2ia = 2.6323E22 * 21--254E-24 = 0.03858 cm' 


-1 

Eqn (3) 

Eroro Eqns. ( 2 ) & (3), we can evolve lilirou^ Egn (l) 


xz 


3 * (0.03858 + 0.54459 * 0.7488) 


0.74674 cm' 
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AS. 2 INCOMPLETE ERROR FUNCTIQS* 


J-i J JU 




n 


j 


- tx 


-xn 




dt = x' 


n-1 


-t 


X 


lH 


dt 


are called incomplete error functions. 
Tney nave series solutions. 


Eq (x) = e”^/x 




In general/ 
oo 


(-X ) 


m 


m*^o 

m^-1 


<n-l-m)ni ‘ 


n n-1 

+ (-1) ' (n.i) j^V' tln(x)-A^) 


= X 

1 n— f — m Jin ’ — — . - 

^ n>- o 
n-1 

where/ f = 0,577216 z A^= 0 / A^= ( ”} 

-»n=l 


In particular ^ 

E^(x) = - ( y +in x) + X- ^ ^ 

2 3 

E2<x) = 1 + x<V+ In x) J2 


32 


The following integral has been used j- 

n- 


J kx n-1 i kx 

e .E^(ax4b),dx = .e % E„ ^ (ax-rt>) - 


h-±' 


i=0 


n-1 bk , 3x4b 

^ E, < p — ) 

TF"*e f 1 

k a 


Ref(31)^ page 355 
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In particular/ 


kx 

e .E^(ax+b).dx 


> 10 


f — E (ax+b) - , _ ■ ax-l-b , 

"" 1 k 1-k/a ' 

E^(x) is well approximated by;- 


E (x) = e“^ 
n n 

1 ^ 

where/ Q (x) =-— -7 — ; — r 

n x+n (x+n; 


A6.3 


(a) Solution ; 1*^ 0, (r) - c^*0. (r) = 


-Ir 

pg(r) - c^* 0^(r) + P % 


In sflierical coordinates 


= -i— * ) 

^ 2 dr ' ^ * dr ' 


where P - Qe 


i ^ 2?_£s. > - c" 

dr dr ® r^ 


Substituting r,^_,(r) = X* # equation reduces to 


(D^-C^) X = - P s 

r 


- rZv 


By using the simple integrating factor method by first 
00 n sidering (D+c) X=Y / solving for Y and later/ 
solving for X, we get the equation 
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-2cr 

X = -4^ + 

-2c 

From appenaix A6.2, 
As X = 


cr r -2cr f 

- P.c J dr.o J — .c3r 
. cl” = - E^( (c-2I ).r ) 


/ * r 


Ae 


-2cr 


Be 


cr 


2c 


+ ?e 


-Jdr. 


— 2cr 


.Ej.Cl'.tofSJ) 


Again using the integral identity of A6.2, 




-2cr 


Be 


cr 


-tr^e 


cr 


(c-S>) 
2c 


■ r ^ 2c ^^1^“ ^,2c 


(c-Z 


c-r 


As 0g(r) 


►- 0 as r — yinf , B = O is essential. 
Thus, substituting back for P, we obtain:- 


— 2 dr jj "N 

0jr) = ^{c-2)*E^(Y)e'^ “ ®l{X)*e""''^ J 

v^ere, X = (c-^),r 

Y = >7(1+ 

(b) Solution of 0 (r) - c^. 0 (r) + Qe“^ = 0 
— '^S s 

Here the approacii is similar. Substituting for 


r * ^(^) = X^ we obtains- 
(P^-c^) X = -Pr.e” 


V 


Eqn (1) 
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The complementary Function of i :iis differential Eqn. is 

-cr cr 

X = Ae + Be 

The particular Integral is obtaineo by choosing a particular 
solution for X = (Dr + E)e”’ ; 


Putting this is eqn. (l), we obtain the coefficients 
1 ^ 2 


D = 




/ E 




-Sr 


—ex’ cx 

Hence, ^ ^ + Q ® ( r - 2 ) 

® ^ F“ c^£^ (c2_j.^)^ 

as 0 (inf) =9, B=0 

O 


• (r) = ~ 

S X 


-cx 






(c^-2?) 


_ ) 
2 > 




TOTAL GAMMA EMISSIONS FROM CX)kE 
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A 7.1 EVALUATION OF MODERAT(m TIME CONSTANT 


Heat balance for a fuel rod channel gives 
leat added) = (heat suppliec to co< Lant) - ( heat carrl.ed away 
flowing coolant ) . 


r) mCp '^ac 
^ dt 


Ah (T ) + ft L (T.-T ) .... Eqnd) 

c cl ac p i ac ^ 


where, Tj_ , T = coolant inlet/average teinprs. in core 

3.0 

= Cladding outer tempr. 

A = coolant flow rate 


jninatiny ^cl in Eqn(l) by substitutiny 

got by Steaay-state 


h. T . + h T 
T - £ af c ac 

■^cl h_ + h 


t transfer between fuel & coolant, indiere h^, are fuel 

1 transfer coefficient and average tempr. given by’'Eqn(7. 17) / 
>btains- 


d T . 
ac 

dt 


A h 




( T -T ) - ~ 
af ac m 


-<T -T.) . . Ectn(2) 

ac i 


e moderator time copsl^ant is the coefficient of its differential 


Eqn (3) 

5 evaluated from Eqri(2) by' putting *^'^ac 

dt ~ ° 

* -J 

tr steady State condition). We obtain :t 


1 



